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Preface 


This  structural  model  investigation  on  the  Montgomery  Point  torque-tube 
gates  was  authorized  by  Headquarters,  U.S.  Army  Corps  of  Engineers 
(HQUSACE),  on  15  October  1993  at  the  joint  request  of  the  U.S.  Army  Engi¬ 
neer  Division,  Southwestern  (USACE-SWD),  and  the  U.S.  Army  Engineer 
District,  Little  Rock  (USACE-SWL).  The  studies  were  conducted  by  person¬ 
nel  of  the  Hydraulics  (HL)  and  Structures  Laboratories  (SL),  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  during  the  period  November 
1993  to  February  1995  under  the  direction  of  Mr.  F.  A.  Herrmann,  Jr., 
Director,  HL;  Mr.  B.  Mather,  Director,  SL;  Mr.  R.  A.  Sager,  Assistant 
Director,  HL;  Mr.  John  Ehrgott,  Assistant  Director,  SL;  and  under  the  gen¬ 
eral  supervision  of  Mr.  G.  A.  Pickering,  Chief,  Hydraulics  Structures  Divi¬ 
sion  (HSD),  HL;  Dr.  J.  P.  Balsara,  Chief,  Geomechanics  and  Explosion 
Effects  Division  (GEED),  SL;  Dr.  R.  Mosher,  Chief,  Structural  Mechanics 
Division  (SMD),  SL;  Mr.  N.  R.  Oswalt,  Chief,  Spillways  and  Channels 
Branch,  HSD;  and  Dr.  R.  L.  Hall,  Chief,  Structural  Analysis  Group,  SS-A, 
SL.  Project  Engineers  for  the  model  studies  were  Mr.  B.  P.  Fletcher,  HSD, 
and  Dr.  L.  A.  de  Bejar,  SS-A.  Mrs.  S.  Garner,  SMD,  assisted  with  the  finite 
element  analysis  of  the  gate.  Instrumentation  support  was  provided  by 
Messrs.  H.  C.  Greer  III,  T.  W.  Warren,  J.  C.  Abies,  and  W.  Guy,  Instru¬ 
mentation  Services  Division,  Information  Technology  Laboratory.  This 
second  report  of  the  investigation  was  written  by  Dr.  de  Bejar. 

During  the  structural  investigation,  Messrs.  L.  Winters  and  T.  Papageorge 
(USACE-SWL-ED-DS)  and  Mr.  R.  Veselka  (USACE-SWD-ED-TS)  visited 
WES  to  discuss  the  structural  models  and  provide  valuable  suggestions. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes,  Gtation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows; 


Multiply 

By 

To  Obtain 

cubic  feet  per  second 

0.02831685 

cubic  metres 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres 

inches 

25.4 

millimetres 

kips 

4.448222 

kilonewtons 

kips  per  foot 

1,355.818 

newton-metres 

kips  per  square  foot 

47.88026 

kilopascals 

kips  per  square  inch 

6.894757 

megapascals 

pounds 

4.448222 

newtons 

pounds  per  cubic  foot 

16.01846 

kilogram-force  per  cubic  metre 

pounds  per  foot 

14.5939 

newtons  per  metre 

pounds  per  square  foot 

47.88026 

pascals 

pounds  per  square  inch 

0.006894757 

megapascals 

slugs 

14.5939 

kilogram-mass 
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Introduction 


Montgomery  Point  Lock  and  Dam  will  be  constructed  near  the  mouth  of 
the  White  River  Entrance  Channel,  which  is  the  passageway  from  the  Missis¬ 
sippi  River  to  the  McClellan-Kerr  Arkansas  River  Navigation  System.  The 
flow  through  the  navigable  pass  spillway  will  be  controlled  by  a  parallel  array 
of  10  torque-tube  gates.  Due  to  complex  hydraulic  loading,  the  gates  are 
subjected  to  stressing  steady-state  flow  and  random  vibrations. 

This  is  the  second  of  a  two-part  report  on  a  research  project  to  assist  the 
engineer  in  designing  the  torque-tube  gates.  The  objective  of  the  first  report 
is  to  provide  the  designer  with  experimentally  based  information  on  the  struc¬ 
tural  response  of  the  gate  system  in  operation.  This  second  report  describes 
several  mathematical  models  to  represent  the  effect  of  the  hydraulic  pressure 
field  on  the  structure  of  a  typical  Montgomery  Point  torque-tube  gate  and  its 
experimental  verification  from  tests  on  a  corresponding  1/15-scale  physical 
model  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES). 


Purpose  and  Scope  of  the  Research  Project  on  the 
Structure 

The  research  project  on  the  structure  satisfies  multiple  purposes: 

a.  Construct  and  instrument  a  rigid  1/ 15-scale  physical  model  of  the  gate 
to  experimentally  determine  (1)  reaction  forces  for  prototype  design  of 
the  support  bearings,  and  (2)  time-histories  of  hydraulic  pressure  at 
select^  discrete  locations  on  both  sides  of  the  gate  skin. 

b.  Determine  a  design  representation  of  the  net  hydrodynamic  pressure 
field  acting  on  the  gate  under  (1)  regular  operating  conditions,  and 
(2)  extreme  conditions  of  discharge  and  pool  elevations,  and  unusual 
operation,  or  both.  These  semiempirical  representations  of  the  pres¬ 
sure  field  on  the  gate  are  based  on  sound  theoretical  principles  and 
adjusted  to  fit  the  results  from  numerous  automatically  controlled  tests. 
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c.  Construct  an  analytical  model  of  the  gate  for  the  finite  element  analysis 
of  the  structure  subjected  to  the  design  hydraulic  pressure  field.  This 
analysis  renders  design  reaction  forces  at  the  support  bearings  to  be 
compared  with  the  experimental  measurements. 

d.  Determine  the  gate  natural  frequencies,  mode  shapes,  and  parameters 
of  energy-dissipation  mechanisms  using  (1)  finite  element  analyses, 
and  (2)  nondestructive  experimental  techniques.  Knowledge  of  the  gate 
eigenproperties  is  necessary  for  effective  assessments  of  structural 
safety  during  the  useful  life  of  the  gate. 

e.  Develop  a  simple  structural  model  for  engineering  analysis  of  the 
dynamic  response  of  the  gate  under  hydraulic  flow.  A  random- 
vibration  model  for  the  stochastic  description  of  the  fluctuations  of  both 
the  hydraulic  pressure  field  and  selected  parameters  of  structural 
response  about  the  respective  steady-state  mean  components  follows  a 
simple  deterministic  representation  of  the  gate  system  subjected  to  the 
hydraulic  forces. 

f.  Evaluate,  by  spectral  analysis,  the  safety  of  a  typical  torque-tube  gate 
in  service  against  fatigue  from  flow-induced  vibrations  over  its  design 
life. 


Organization  of  Report  2 

Chapters  in  this  report  are  organized  in  a  prerequisite-development 
sequence.  Thus,  the  first  step  is  to  develop  an  analytical  representation  for 
the  hydraulic  forces  on  the  gate.  In  Chapter  2,  the  steady-state  mean  pressure 
field  is  derived  from  a  semiempirical  application  of  the  theorem  of  moment- 
of-momentum  in  fluid  dynamics.  This  theoretical  prediction  tells  the  analyst 
what  to  expect  as  input  during  the  tests  (Photos  1  trough  3  exhibit  the  series 
of  three  gates  in  the  experimental  flume  at  WES)  and  suggests  the  general 
trend  for  a  proposed  design  pressure  envelope  from  the  results  of  Phase  1  of 
the  experimental  program,  which  involves  the  regular  operating  conditions  of 
the  gate. 

The  experimentation  follows  together  with  finite  element  analyses.  The 
results  of  these  experimental  analyses  (Chapter  3  and  the  second  part  of  Chap¬ 
ter  4)  support  the  final  design  of  the  gate  system.  Chapter  3  describes  the 
details  of  the  experimental  setup  and  the  results  of  finite  element  analyses  for 
both  the  physical  model  (a  brass  structure  standing  on  three  supports)  and  the 
prototype  gate  system  (made  of  structural  steel  standing  on  five  supports). 
Reaction  forces  for  prototype  design  of  the  support  bearings  are  obtained 
through  finite  element  analysis  of  the  prototype  subjected  to  the  scaled  design 
pressure  field  on  the  physical  model.  The  reaction  forces  from  the  finite  ele¬ 
ment  analysis  provide  the  input  for  the  final  design  of  bearings. 
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In  Chapter  4,  two  separate  models  use  the  experimental  data  from  Phase  1 
of  the  experimental  program  to  estimate  extreme-value  wide-hand  hydrody¬ 
namic  pressure  fluctuations  during  the  useful  life  of  the  gates,  with  excellent 
agreement  between  the  two  results.  The  chapter  ends  with  a  description  of 
Phase  2  of  the  experimental  project,  which  generally  considers  very  large 
hydraulic  pressure  fields  on  the  gate  due  to  extreme  conditions  of  discharge 
and  pool  elevations,  and  unusual  operation,  or  both. 

In  Chapter  5,  natural  frequencies  and  mode  shapes  obtained  from  different 
experimental  techniques  (Twang  tests,  forced-vibration  tests,  and  Laser  vibro- 
meter  tests)  applied  to  the  data  are  successfully  compared  with  the  correspond¬ 
ing  analytical  results  from  finite  element  analyses  on  the  grid  for  the 
small-scale  physical  model.  Knowledge  of  these  eigenproperties  is  necessary 
for  effective  assessments  of  structural  safety  during  the  useful  life  of  the  gate. 

Chapter  6  discusses  the  approximate  evaluation  of  the  natural  frequency  of 
free  vibration  of  the  gate  and  the  identification  of  the  equivalent  viscous 
damping  ratio  for  the  physical  model  (both  dry  and  upstream  submerged). 

Two  ramifications  of  an  engineering  model  for  the  dynamic  response  of  the 
torque-tube  gate  system  are  presented.  The  first  section  describes  the  struc¬ 
ture  dynamics  model  from  an  entirely  deterministic  point  of  view.  The  gov¬ 
erning  equation  of  motion  for  the  gate  system  is  formulated  and  solutions  are 
derived  for  the  mean  of  selected  parameters  of  response.  The  second  section 
describes  a  random  vibration  model  for  the  stochastic  description  of  the  fluctu¬ 
ations  of  both  the  hydraulic  pressure  field  and  selected  parameters  of  struc¬ 
tural  response  about  their  steady-state  mean  components. 

The  results  of  the  simulations  presented  in  Chapters  4  through  6  suggest 
further  investigations  on  extreme-value  prediction,  gate  parameter  identifica¬ 
tion,  and  practical  engineering  estimates  of  structural  response. 

Chapter  7  outlines  the  conclusions  of  the  study  and  formulates  recommen¬ 
dations  for  design. 

The  report  ends  with  an  elemental  spectral  analysis  to  assess  the  fatigue 
performance  of  the  gate  system  in  service  over  the  long  run.  A  set  of  selected 
power  spectral  density  functions  for  pressure-cell  records  in  the  controlling 
experiments  illustrate  some  of  the  salient  points  in  the  discussion. 


Chapter  1  Introduction 


3 


Chapter  1  Introduction 


Photo  1.  Montgomery  Point  Torque-Tube  Gate  Model.  Aperture  in  experimental  flume  corresponding  to  central  (instru¬ 
mented)  gate  in  fully  open  position 


experimental  flume  in  fully  open  position 
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Photo  3.  Montgomery  Point  Torque-Tube  Gate 
flume  in  fully  closed  position 


2  Mathematical  Model  for 
Hydraulic  Pressure  Field 


The  study  begins  by  introducing  a  semiempirical  model  for  the  hydraulic 
pressure  field  on  the  gate.  The  objective  of  this  representation  is  to  predict 
the  hydraulic  forces  on  the  gate  in  service.  The  model  will  allow  expedient 
estimation  of  forces  on  the  prototype,  for  engineering  purposes,  as  well  as  on 
other  small-scale  torque-tube  gates  in  future  experimental  projects. 

Water  is  assumed  an  incompressible,  inviscid,  and  homogeneous  fluid; 
and  its  flow  is  taken  as  irrotational.  The  critical  hydraulic  effect  on  the  struc¬ 
tural  model  is  considered  to  occur  when  the  upstream  flow  at  the  maximum 
operational  headwater  pool  runs  over  the  torque-tube  gate  in  a  fully  raised 
position.  This  critical  condition  is  simplified  in  the  model  by  neglecting  the 
effect  of  the  tailwater  pressure. 

This  idealization  leans  on  experimental  evidence.  During  the  experiments 
with  the  physical  model  at  WES,  the  maximum  reaction  response  for  regular 
operation  was  obtained  when  a  unit  discharge  of  water  of  50  cfs  runs  over  the 
system  of  fiilly  raised  gates  (70-deg  elevation  (el))  with  the  tailwater  pool 
maintained  at  95  ft  (prototype  level)  and  the  headwater  pool  growing  mono- 
tonically  from  about  114  to  119  ft  (prototype  level). 

Plate  1  shows  a  lateral  view  of  the  physical  model  of  the  gate  in  closed 
position  (70-deg  el)  with  a  ventilated  nappe  and  prototype-equivalent  95-ft 
tailwater  pool  (the  model  described  in  this  section  will  be  developed  in  the 
context  of  prototype  dimensions).  The  thickness  of  the  nappe  root  (water 
depth  over  the  gate  ridge)  is  carefiilly  measured  to  scale  in  this  picture  as  Ay  = 
2.9  ft.  Similarly,  the  nappe  thickness  at  the  level  of  the  midheight  of  the 
approaching  flow  is  (approximately)  determined  experimentally  as  d'  =  1.6  ft. 

Plate  2  shows  a  plan  view  of  the  geotechnical  borings  at  the  site 
(U.S.  Army  Engineer  District,  Little  Rock  (USACE-SWL)  1993).  The  corre¬ 
sponding  survey  determines  the  natural  elevation  of  the  White  River  bed  at 
selected  locations  on  the  upstream  side  of  the  projected  navigational  pass. 
Standard  Penetration  Test  boring  station  No.  83  at  200  ft  from  the  gates  is 
particularly  noticed  at  a  level  of  88.1  ft  in  Plate  3.  This  location  is  also 
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schematically  located  in  Plate  4  exhibiting  a  longitudinal  cross  section  of  the 
prototype  system. 

Assuming  a  parabolic  upper  surface  for  the  nappe,  its  equation  is  easily 
determined  from  the  condition  of  fitting  the  thickness  of  the  nappe  root,  i.e., 
h  =  hpdX  the  location  of  the  gate  ridge.  The  following  simple  expression  is 
obtained: 


— T 

1, 15.03  j  ’ 


(1) 


where  both  x  and  y  are  in  feet. 

Using  this  expression  and  the  geometry  of  the  circle  with  center  at  C  in 
Plate  4,  toe  following  equation  for  toe  general  water  depth  over  toe  gate  is 
derived: 


h  =  H  -  asin0  -  2.346  (1  - 1.0643  cos0)^ 


(2) 


where 

H  =  27.4  ft  is  toe  level  difference  between  points  O  and  C 

a  =  24.5  ft  is  toe  radius  of  toe  circle 

6  £  (20  deg,  72.75  deg)  =  elevation  of  a  general  point  K  on  toe 
gate  skin 

Recognizing  that  the  streamlines  going  through  toe  sections  at  G  and  G'  of 
the  control  volume  between  them  in  Plate  4  are  locally  approximately  parallel 
to  each  other  and  that  toe  gauge  pressure  at  G'  is  practically  zero,  toe  follow¬ 
ing  Bernoulli  equation  is  formulated  for  streamline  GG'  (Prandtl  and  Tietjens 
1957a, b.  Shames  1962): 


(3) 

'2-g  y  2g  ’ 


where 


Pq  =  gauge  pressure  at  point  G 

Vq  and  Vq,  =  assumed  uniform  flow  velocities  at  points  G  and  G', 
respectively 
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7  =  unit  weight  of  water 
g  =  acceleration  of  gravity 

Invoking  the  principle  of  continuity  between  the  cross  sections  at  G  and 
G'  (Prandtl  and  Tietjens  1957a, b,  Shames  1962),  one  may  solve  for  the  veloc¬ 
ity  Vq  to  get 


^a  = 


gd 


(4a) 


or,  evaluating  this  expression  (using  d  =  30.9  ft,  the  depth  of  the  approaching 
flow). 


Vq  =  1.64  fps 


(4b) 


This  value  of  the  approaching  flow  velocity  compares  very  well  with  the 
value  obtained  from  the  design  approaching  unit  flow  =  50  cfs,  i.e.. 


1.62  fps  . 


(4c) 


The  average  value  Vq  =  1.63  fps  will  be  used  in  the  development  of  the 
model. 

Notice  that  although  streamlines  at  cross  section  OA  at  the  base  of  the 
gate  are  not  strictly  parallel  to  each  other,  the  Bernoulli  equation  for  line  GC?’' 
provides  an  approximate  estimate  for  the  total  pressure  head  at  G''  (Olson 
1961,  Vennard  and  Street  1975),  giving 


=  9.43  ft  , 
Y 


which  is  close  to  the  hydrostatic  pressure  head  at  the  same  location  (9.5  ft). 

Likewise,  since  the  streamlines  are  practically  parallel  to  each  other  at 
cross  sections  G'  and  G",  respectively,  i.e.,  at  the  ends  of  the  nappe  portion 
being  considered,  the  Bernoulli  equation  for  line  G'G"  provides  an  accurate 
estimate  for  the  hydrodynamic  pressure  head  at  G",  giving 
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=  -3.573  ft  . 
Y 


If  the  total  pressure  head  at  a  point  is  envisioned  as  the  summation  of  hy¬ 
drostatic  and  hydrodynamic  components  (Prandtl  and  Tietjens  1957a, b),  i.e., 


£  =  P.  ^  R 
Y  Y  Y  ’ 

and  the  hydrodynamic  pressure  head  is  assumed  to  vary  linearly  from  G''  to 


from  ^  =  -0.07  ft  to  ^  =  -3.573  ft  . 
Y  Y 


recognizing  that  the  hydrostatic  component  at  a  general  point  G'"  is  one-half 
of  the  corresponding  water  depth,  the  following  expression  for  the  total  pres¬ 
sure  head  at  (?'"  is  obtained,  as  a  function  of  the  point  elevation  (Plate  5): 


=  14.95816-12.25sin0-1.173(l-l.O643cos0)2-O.O664O80  .  (5) 
Y 

Therefore,  the  general  energy-balance  equation  for  line  GG"''\n  Plate  4 
provides  an  approximate  means  to  determine  the  water  depth,  h,  at  any  point 
G'"  out  directly  over  the  gate  AB,  i.e.. 


/f-  a  sin  0  +  — 
2^ 


(df  h  m=o 
UiJ  2  Y 


(6) 


where  p(6)  is  given  by  Equation  5. 

Simultaneous  solution  of  Equations  5  and  6  allows  the  numerical  evalua¬ 
tion  ofh,  resulting  in  a  negligible  variation  from  the  parabolic  profile  origi¬ 
nally  assumed. 

Consider  now  the  idealized  unit-thickness  control  volume  D 'DA*'' ABB' 
A"A'"D'  in  Plate  6  and  its  natural  partition  into  two  components  by  the  ver¬ 
tical  plane  A" 'A*''.  The  approaching  flow  may  be  assumed  to  enter  and  exit 
the  volume  with  uniform  velocities  U  and  u,  respectively.  The  forces  acting 
on  the  boundary  surface  are  schematically  represented  in  the  diagram.  The 
weights  of  water  in  the  various  parts  of  the  control  volume  are  designated  by 
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W;  i  =  1,2,3.  The  vertical  and  horizontal  components  of  the  hydraulic  force 
exerted  by  the  flow  on  the  gate  are  designated  by  and  respectively. 

The  plate  shows  the  components  of  the  corresponding  reactive  force,  i.e.,  the 
force  exerted  by  the  gate  on  the  fluid,  after  sliding  along  its  line  of  action 
from  the  point  of  application  Q'  to  point  Q  on  the  gate  skin  surface  (Plate  7). 
These  forces  have  lever  arms  r,  and  r2  with  respect  to  the  pivoting  point  A, 
respectively. 

Assuming  that  the  flow  through  AA"  is  approximately  horizontal,  the  ver¬ 
tical  equilibrium  of  control  volume  ABB' A' 'A  may  be  expressed  as  (Prandtl 
and  Tietjens  1957a,b,  Shames  1962) 

7?,  -  1^2  -  1^3  =  0  (7) 


In  the  absence  of  horizontal  body  forces,  the  equation  of  conservation  of 
linear  momentum  applied  to  the  whole  control  volume  in  Plate  6  may  be 
expressed  as 

=  p  p 


where  and  P2  represent  the  resultant  forces  of  pressures  acting  on  sides 
DD'  and  AA*''  of  the  control  volume,  respectively;  Aj^dxiAAg  are  the  cross- 
sectional  areas  of  the  control  volume  at  locations  D  and  B,  respectively;  U  and 
u  are  the  corresponding  flow  velocities,  assumed  to  have  a  direction  approxi¬ 
mately  normal  to  the  respective  cross-sectional  areas;  and  p  is  the  mass  den¬ 
sity  of  water. 

The  theorem  of  moment-of-momentum  about  point  A,  applied  to  the 
whole  control  volume,  produces  the  following  equation: 


H 

—  -  a  sinflr 


+  (COS0O  -  cos0^)  +  +  ^2/2 


-  7^r2  =  -|  y  -  asin^o  ^if-pAjj 


hf 

a{s\ndf  -  sin^o)  +  — 


u^pAg  , 


(9) 


where 

77  and  Ay  =  water  depths  at  sides  DD'  and  BB'  of  the  same  control 
volume,  respectively 
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Bq  and  Bj-  =  elevation  angles  of  points  A  and  B,  which  are  the  respective 
ends  of  the  circular  gate  AB  in  closed  position  (70-deg 
elevation  with  respect  to  the  horizontal  plane) 

a  =  circle  radius  (Plate  4). 

Since  point  Q  is  on  the  circle,  the  following  relationship  is  satisfied; 

(acos^o  “  +  (asin^o  +  =  a^  .  (iO) 


Next  consider  the  generating  radius  vector  with  center  at  C,  magnitude  a, 
and  argument  B  shown  in  Plate  8.  The  shaded  area  is  given  by 


jcos^o  JtsinSf 

f  f  dydx 


3COS0.  J~2  ~Z2 


= 


1 

sin0/cos0o-cos0y)  +  —  (sin20^  -  sin20o)  -  — - — 

4  12 


(11a) 


/j 


and  its  centroidal  point  has  abscissa 


;>cos0Q 


J  asinflo 


2A, 


^(sin^Bf-  sin^0o)  “  a^sin^^o  (sin0y^  -  sin^o) 


(lib) 


Therefore,  the  weight  is  obtained  from 


W,  =  yA^ 


(11c) 


where 

7  =  pg  is  the  unit  weight  of  water, 
and  the  location  of  its  line  of  action  can  be  calculated  from 
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Aj  =  acos^o  "  ^ 


(lid) 


On  the  other  hand,  the  resultant  forces  P,-,  i  =  1,2,  can  be  estimated  by 
the  expressions  (Prandtl  and  Tietjens  1957a, b) 


p  - 

*  ~T~ 


(12a) 


and 


1  -  — sin0 
H 


(12b) 


corresponding  to  the  hydrostatic  pressure  distributions  on  DD'  and  A  A*'', 
respectively. 

Evaluating  for  the  prototype  gate  dimensions,  a  =  24.5  ft  and  //  = 

27.4  ft,  with  Bq  =  20  deg.  By  =  72.75  deg  and  7  =  62.4  pcf.  Equations  11a 
and  11c  render 


=  80.91  ft^  (13a) 

and 

W'  =  5.049  k  ,  (13b) 


respectively;  whereas  Equations  1  lb  and  1  Id  render 


X  =  18.66  ft 


(13c) 


and 


=  4.36  ft 


(13d) 


respectively. 

Likewise,  Equations  12a  and  12b  yield 
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(14a) 


Pj  =  23.424  k 


and 

Pj  =  12.136  k  , 


(14b) 


respectively. 

Now  the  reactive  force  Pj  may  be  computed  from  Equation  7  as 

Pi  = 

=  y  ah (cos9q  -  cos6^  + 

=  2.852  P+  5.049;!:=  7.901  . 

The  experimental  water  depth  at  the  root  of  the  nappe  was  previously 
obtained  as 


hf  =  2.9  ft  .  (16) 

and  the  flow  velocities  U  and  u  may  be  estimated  from  direct  application  of 
the  equation  of  continuity,  i.e., 


£/  = 


=  1.63  fps 


30.9 

27.4 


1.84  fps 


and 


rn  H 
u-U\  — 


hf 

-  1.84  fps 


^  j  =  17.39  fps 


(17a) 


(17b) 
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Putting  Equations  14a,  14b,  16,  17a,  and  17b  into  Equation  8,  the  reactive 
force  P2  may  also  be  computed  as 
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(18) 


g 

=  9.768  k  . 

Plate  9  shows  the  free-body  diagram  of  the  idealized  gate  subjected  to  the 
hydraulic  forces  Ri  and  R2  (per  unit  width  of  gate)  and  the  corresponding 
internal  reactions  at  the  base  of  the  cantilever,  where  Nq,  Vq,  and  Mq  repre¬ 
sent  the  axial  load,  shear  force,  and  bending  moment  per  unit  width  of  gate. 
The  equilibrium  formulation  leads  to  the  following  force  transformation 
(Boresi  and  Sidebottom  1985,  Przemieniecki  1985): 


'^0 

►  — 

1 

COS0O 

sin0Q 

-^1 

,^0. 

-sin^o 

cos^o 

(19a) 


Equation  19a  allows  the  direct  evaluation  of  the  shear  force: 


Vq  =  i?jsin0o  +  R^cosBq  ^ 

=  11.881  vlr 

The  result  in  Equation  19b  provides  one  of  the  boundary  conditions  for 
the  boundary-value  problem  to  be  formulated  next.  The  other  boundary  con¬ 
dition  is  given  by  the  fact  that  the  curved  cantilever  has  a  free  end  at  B,  i.e.. 


K^o)  =  lo  =  11-881^ 
He,)  =  0 


(20) 


Now,  if  the  analysis  leading  to  Equations  15  and  18  for  the  hydraulic 
forces  is  generalized  for  an  end  elevation  6,  instead  of  ^  (Plate  4),  the  follow¬ 
ing  expressions  are  obtained 


R^iO)  =  Y 


a/i(cos0o  “  COS0) 


+ 


cos0o(sin6 


-  sin^o) 


-  i  (sin20  -  sin20o)  -  Ue  -  6^) 


(21a) 


Chapter  2  Mathematical  Model  for  Hydraulic  Pressure  Field 


15 


where 

h  =  water  depth  of  the  nappe  root  at  the  crest  of  the  generalized 
cantilever  AK. 

After  some  algebraic  manipulations,  the  first  two  derivatives  of  each  of  these 
expressions  may  be  obtained  as 


- 1  =  va  COS0 

de 


h\  a  T7^(cos0o-cos0) 

—  tan0  +  - 

a)  7?-l 


(22a) 


^-^2  ,  T  <  Tf  COS^d 

dO-  [77^  -  1  [  17^  -  1 


3t?~cos^0  I 

(7?2  -  1)3  J 


where 


V 


2 


(22d) 


(22e) 


Next,  consider  the  free-body  diagram  of  the  differential  element  of  gate  in 
Plate  10.  The  equilibrium  of  forces  in  the  radial  and  tangential  directions, 
respectively,  leads  to  the  following  relations: 
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y?/sin0  +  ^2^cos0  +N=  — 

Ciu 

and 

RJsvnO  -  Ji,'cos6  +  V  +  -^  =  0  ,  (23b) 

dd 

where  the  prime  symbol  represents  derivative  with  respect  to  6. 

Eliminating  dNIdd  from  these  equations,  one  obtains  the  following  govern¬ 
ing  differential  equation: 


IR'cosO  -  IR'sLne  +  ;?/'sin0  +  R," cosO  =  ^  +  V  ,  (23c) 

which  together  with  the  boundary  conditions  (Equation  20)  constitute  a 
boundary-value  problem  for  the  shear  force  function  V{B)  (Kreyszig  1988). 

To  facilitate  the  numerical  solution  of  Equation  23c,  the  boundary-value 
problem  will  be  transformed  into  an  initial-condition  problem  (MATHSOFT 
1994).  Since  V(0)  is  already  known,  only  V'(0)  is  missing  to  complete  the 
formulation.  The  problem  is  entirely  solved  with  MATHCAD  5.0  (see  Pro¬ 
gram  1:  “Theoretical  Estimate  of  Pressure  Field,”  at  the  end  of  this  chapter). 
Initially,  a  guess  value  for  V'(0)  is  provided.  Subsequently,  the  program 
'sweeps'  the  6-domain  numerically,  evaluating  the  function  V(6)  until  the  value 
V(6y)  at  the  other  end  of  the  domain  is  obtained.  If  this  value  coincides  with 
that  one  stipulated  by  the  original  boundary-value  problem  [V(6^)  =  0],  within 
a  previously  established  tolerance,  then  the  starting  initial  condition  V'(0)  is 
correct  and  the  initial-condition  problem  is  fully  formulated.  Otherwise,  a 
new  trial  value  is  automatically  advanced,  and  the  process  is  repeated  iter¬ 
atively  until  convergence  is  reached. 

Once  the  value  of  V'(0)  is  obtained,  the  governing  differential  Equa¬ 
tion  23c  is  solved  using  the  Fourth-Order  Runge-Kutta  method  for  initial- 
condition  problems  (Kreyszig  1988).  The  solution  function  =  V(6) 

and  its  derivative  Z'^“^  =  dVIdd  are  plotted  within  the  hardcopy  of  the 
workspace  for  the  submitted  task.  Notice  that  the  hydraulic  pressure  field  is 
given  by 


m  = 


1 

a  dO 


(24) 


and  appears  plotted  and  listed  both  for  the  prototype  gate  and  for  the  model  at 
the  end  of  the  workspace,  with  the  objective  of  direct  comparison  with  the 
experimental  pressure  field  on  the  small-scale  physical  model. 


Chapter  2  Mathematical  Model  for  Hydraulic  Pressure  Field 


17 


Plate  11,  generated  with  the  program  MATLAB  4.0  (MATHWORKS 
1993),  shows  the  minimum  and  maximum  values  of  the  experimental  pressure 
acting  on  the  upstream  side  of  the  gate  during  the  controlling  100-s  test 
(T1080).  The  pressure  field  is  described  along  three  different  gate  arcs:  the 
left  edge,  the  middle  arc,  and  the  right  edge  (Plate  12).  Notice  that  the  curves 
represent  polynomial  fits  on  three  data  points,  which  explains  the  intersection 
between  the  upper  and  the  lower  bound  curves  for  the  middle  arc  toward  the 
gate  base.  These  minimum  and  maximum  curves  delimit  narrow  bands,  sug¬ 
gesting  very  little  fluctuation  of  the  pressure  at  a  given  location  during  the 
test.  The  ftieoretical  estimation  of  the  pressure  field  is  an  excellent  upper 
bound  over  most  of  the  domain,  except  on  a  small  region  toward  the  gate 
base,  where  the  estimation  is  slightly  unconservative. 

The  assessment  of  the  hydraulic  moment  Mq  requires  the  previous  evalua¬ 
tion  of  the  lever  arms  Tj  and  r2.  These  lever  arms  are  obtained  from  the 
simultaneous  solution  of  Equations  9  and  10.  Notice  that  evaluation  of  the 
segment  y2  is  required  to  define  this  system  of  equations.  Segment  y2  may  be 
computed  from 

72  =  asin^o  -  7  , 


where 


_L 

P,  6 


2 

- 

1  -  — sin^o 

sin^o 

H  "J 

.(25b) 


This  system  of  equations  is  solved  numerically  using  the  program  MATH- 
CAD  5.0  (MATHSOFT  1994)  (see  Program  2:  “Location  of  Resultant 
Hydraulic  Load,”  at  the  end  of  this  chapter),  yielding  the  values 


r,  =  3.434  ft 


(26a) 


and 

^2  =  6.336  ft 


(26b) 


Consequently,  the  hydraulic  moment  per  unit  width  of  gate,  Mq,  may  be 
obtained  as 


+  P2^2 

=  89A-ft  . 


(27) 
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The  theoretical  total  hydraulic  forces  on  a  prototype  gate  (total  width  = 

30  ft)  appear  listed  in  Table  1.  The  experimental  counterpart  values 
(USAEWES  1994a)  were  obtained  from  data  on  ventilated-nappe  tests  T1087 
(sampling  rate:  100  sps;  duration:  500  s)  and  T1 109  (sampling  rate: 

500  sps;  duration:  200  s).  The  experimental  forces  are  scaled  up  from  the 
model  by  multiplying  by  a  factor  =  15^.  The  hydraulic  moment  is  scaled 
up  by  a  factor  X^'^  =  15'^.  In  general,  there  is  good  agreement  because  the 
experimentally  measured  moment  is  really  bending  moment  at  the  root  of  the 
operating  arms  not  at  the  torque-tube  center,  and  the  theoretical  vertical  com¬ 
ponent  of  reactive  force  (Ri)  does  not  include  the  effect  of  the  hydrostatic 
pressure  due  to  the  tailwater  acting  on  the  downstream  side  of  the  gate 
(Plate  1). 


Table  1 

Theoretical  Prediction  versus  Experimental  Measurements 

Hydraulic  Forces  on  70  deg  gate  (Tailwater  pool:  95  ft)  | 

Vertical  Force,  k 

Horizontal  Force,  k 

Hydraulic  Moment,  k-ft 

T1087 

172 

313 

2,466 

T1109 

172 

312 

2,350 

i^Bmi 

.237 

'293. 

.2,670 

Average 

Difference 

+  38% 

-6% 

+  11% 
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Plate  1 .  Montgomery  Point  Torque-Tube  Gate.  Lateral  view  of  physical  model  in  closed  position  with  a  venti¬ 
lated  nappe  {T1 109) 
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Plate  2.  Montgomery  Point  Lock  and  Dam.  Plan  of  borings  at  the  site  (Design  Memorandum  No.  1,  Plate  31/1) 
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Plate  3.  Montgomery  Point  Lock  and  Dam.  Boring  profile  C  (Design  Memorandum  No.  1,  Plate  31/3) 


Plate  4.  Montgomery  Point  Torque-Tube  Gate.  Longitudinal  cross  section  of  prototype 
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Plate  6.  Montgomery  Point  Gate  System.  Idealized  control  volume 


.  Shaded  area  subtended  by 
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counterpart  (T1080) 
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USAE-Waterways  Experiment  Station,  Structures  Laboratory 
Vicksburg,  Mississippi 


Montgomery  Point  Gate  Model 
Theoretical  Estimate  of  Pressure  Field 


Dr.  Luis  A.  de  Bejar 
November  9, 1994 


gamma  =62.4 


Unit  weight  of  water 


a  :=24.5 


Radius 


U  :=1.84 


Approach  VELOCITY  in  fps 


g:=32.2 


Gravity 


H:=27.4 


Channel  depth  upstream  from  gate 


:=  19.021 


Initial  value 


6:=  20..  73 


Elevation  Range 


h.  27.4-  24.5-sin —-e  -  2.346-  1  -  1. 0643  -  cos  U^- 
180  180 


Water  depth 


U^\  /H\^ 


S0O  sin  20- 


C0O  cos  20 


theta^  ,^-0- 


s0  :=sin( theta) 


c0  cos( theta) 


oi  2  /h\  C0O-C0 

Rlpr  =  gamma-a  ■  -  •s0-h  c0- - 

1_± 


R2pr  :=gamma  a*  h- 


s20:=2-(s0c0) 


c20  ;=(c0  c0-  S0  S0) 


Program  1.  "Theoretical  Estimate  of  Pressure  Field"  (Sheet  1  of  5) 
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rt  :=((Rlprpr-  2  R2pr)  s0-t-  (R2prpr+  2  Rlpr)  c0) 


Vp:  =  -31500  Guess  value  for  V'(0) 

/- 11881 

load(xl  ,v)  :  = 

\  '"o 


D(ifl,V):=  *  score(x2,V)  :=  Vp 

S  :=sbval(v,  1 , 54, D, load, score) 

S  =-2.769- 10“ 

jum\ 

\-  27690  / 

I 

DD(rt,y)  :  = 

\-yo+>« 

Z  :=rkfixed(y, 0.3491 , 1.26973 ,54, DD) 
i  :=  1..54 


Program  1 .  (Sheet  2  of  5) 
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Program  1 .  {Sheet  3  of  5) 
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Program  1 .  (Sheet  5  of  5) 
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USAE-Waterways  Experiment  Station,  Structures  Laboratory 
Vicksburg,  Mississippi 


Program  2.  "Location  of  Resultant  Hydraulic  Load"  (Continued) 
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Program  2.  (Concluded) 


3  Experimental  and  Finite 
Element  Analyses 


The  objective  of  the  experimental  program  is  to  provide  a  realistic  assess¬ 
ment  of  the  hydraulic  forces  on  the  gate.  Complemented  by  theoretical  mod¬ 
els,  the  experimental  analysis  renders  an  effective  pressure  field  to  be  used  as 
input  in  the  finite  element  analysis  of  the  gate.  This  finite  element  analysis 
provides  the  designer  with  the  response  stress  field  within  components  of  the 
gate,  the  reactions  at  the  support  bearings,  and  the  force  in  the  actuator 
cylinder. 

The  1/15-scale  physical  model  of  the  gate  was  made  of  bronze  for  con¬ 
venience  of  construction.  Since  bronze  has  roughly  one-half  the  modulus  of 
elasticity  of  structural  steel  (American  Institute  of  Steel  Construction  (AISC) 
1986,  Flinn  and  Trojan  1990),  plates  are  easier  to  cut  and  assemble,  thereby 
facilitating  the  overall  construction  process.  Another  advantage  of  bronze  as 
the  construction  material  for  the  model  is  to  have  approximately  the  same 
mass  density  as  structural  steel,  thereby  making  the  rigid  gate  model  close  to  a 
true  Froude’s  model  in  free-fall  experiments  (Novak  and  Cabelka  1981). 

In  the  early  stages  of  the  experimental  program,  due  to  imperfections 
inherent  to  the  construction  process  at  the  shop  and  due  to  warping  while 
welding  some  component  plates  into  the  gate,  the  torque  tube  turned  out  to  be 
slightly  misaligned,  introducing  in  this  way  asynunetry  in  the  final  configura¬ 
tion  of  the  model  at  the  site  of  the  experiments.  To  compound  the  situation, 
it  was  not  possible  to  achieve,  within  a  tolerance  satisfactory  for  experimental 
analysis,  horizontal  alignment  of  the  seven  supports  originally  planned  in  the 
prototype  design,  thereby  introducing  further  asymmetry  in  the  experimental 
records  for  structural  parameters,  in  particular,  for  support  reactions. 

The  effects  of  the  asynunetric  configuration  of  the  model  became  evident 
during  the  preliminary  calibrational  experiments  on  the  gate,  thus  becoming  a 
source  of  concern.  Since  the  hydraulic  forces  on  the  gate  are  expected  to  be 
approximately  symmetrical,  the  corresponding  experimental  response  records 
for  selected  structural  parameters  should  similarly  reflect  symmetry.  The 
problem  was  finally  resolved,  to  engineering  satisfaction,  by  reconstructing 
the  gate  with  extraordinary  attention  to  the  effects  of  uneven  cooling  of 
welded  parts  and  by  having  it  standing  on  only  three  supports.  The  multiple 
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difficulties  found  in  the  construction  of  a  balanced  gate  model  standing  on 
supports  aligned  on  the  horizontal  plane  led  to  a  recommendation  from  the 
WES  engineering  team  to  reduce  the  number  of  supports  in  the  final  prototype 
gate  to  five  (the  construction  process  in  larger  scale  is  more  controllable). 

This  recommendation  was  readily  accepted  by  the  sponsor  (USAEWES 
1994a). 

Consequently,  the  experimental  program  on  the  gate  proceeded  on  a 
model  standing  on  three  supports,  and  the  finite  element  analysis  was  effected 
on  two  analytical  models:  (a)  one  for  the  physical  model  on  three  supports, 
and  (b)  another  for  the  prototype  gate  standing  on  five  supports.  Plate  13 
shows  a  plan  view  of  the  physical  model  and  the  instrumented  bearing  blocks 
2,  4,  and  6,  according  to  the  original  notation  for  seven  supports. 

Plate  14  shows  a  schematic  representation  of  a  typical  load  cell  and  the 
corresponding  bridge  wiring.  The  load  cells  consist  of  two  structural  alumi¬ 
num  tubes  (A1  6061-T6,  std  <f»3/8",  L  =  3/4  in.)  per  support  block  (one  tube 
by  the  upstream  face  of  the  model  and  the  other  tube  by  the  downstream 
face),  bolted  at  both  ends,  i.e.,  to  the  bearing  housing  and  to  the  base,  respec¬ 
tively.  Each  tube  was  instrumented  with  four  strain  gauges,  one  per  quadrant, 
as  indicated  in  Plate  14.  The  0-/180-deg  strain  gauges  were  calibrated  to 
sense  the  normal  load  on  the  tube,  upon  its  axial  distortion,  whereas  the  90-/ 
270-deg  strain  gauges  were  calibrated  to  sense  the  shear  force  at  the  ends  of 
the  tube,  upon  its  flexural  deformation  (lateral  displacement  of  the  bearing 
housing  with  respect  to  the  base).  Individual  tubes  tested  satisfactorily  against 
inelastic  lateral  instability,  even  after  reducing  the  number  of  actual  support 
blocks  to  only  three  (The  Aluminum  Association  1981,  1982).  The  length  of 
the  tubes  was  engineered  to  fit  the  expected  deformations  during  the  experi¬ 
ments  to  the  sensitivity  of  the  strain  gauges.  The  performance  of  the  load 
cells  during  the  tests  was  considered  excellent. 

The  model  gate  was  also  instrumented  with  accelerometers  (PCB  302A)  to 
identify  its  natural  structural  properties  from  the  forced  experimentation  (fre¬ 
quencies  of  natural  vibration  and  corresponding  modal  shapes).  Plate  15 
shows  the  location  of  the  accelerometers  on  the  model.  The  corner  acceler¬ 
ometers  1  and  6  are  triaxial  in  nature  to  provide  information  on  the  lateral- 
torsional  modes.  One  of  the  triaxial  accelerometer  axes  is  normal  to  the  gate 
skin,  while  the  other  two  axes  are  tangential  to  it.  One  of  the  tangential  axes 
is  aligned  with  the  gate  ridge.  Accelerometers  2  through  5  and  7  through  8 
are  uniaxial,  sensing  motion  in  the  direction  normal  to  the  gate  skin. 

One  of  the  most  important  objectives  of  the  project  is  the  experimental 
determination  of  the  hydraulic  pressure  field  on  the  gate  in  the  controlling 
operational  configuration.  To  this  end,  the  model  gate  was  instrumented  with 
nine  pressure  cells  (thermally  insulated  Kulite  XTM-190-10  psia)  on  each  of 
its  faces,  to  provide  information  on  the  spatial  variation  of  the  total  (hydro¬ 
static  plus  hydrodynamic)  pressure  field  on  the  gate.  Plates  16  and  17  show 
the  location  of  the  pressure  cells  on  the  upstream  and  downstream  sides  of  the 
model,  respectively.  Notice  that  the  general  triburary  area  for  a  given 
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pressure  cell  number  on  the  upstream  face  corresponds  with  the  same  pressure 
cell  number  on  the  downstream  side  of  the  gate.  For  example,  pressure 
cell  1  on  the  upstream  side  of  the  gate  senses  hydraulic  pressure  on  a  gate 
area  whose  back  experiences  hydraulic  pressure  sensed  by  pressure  cell  1  on 
the  downstream  side  of  the  gate. 

The  experimental  records  for  each  test  typically  consist  of  the  time-history 
of  response  for  the  parameter  sensed  by  each  transducer  on  the  model,  as  well 
as  the  maximum,  average,  and  minimum  readings  from  the  record.  Table  2 
describes  the  average  values  of  the  experimental  hydraulic  pressure  field  on 
the  upstream  side  of  the  gate  for  a  representative  100-s  test  and  is  introduced 
at  this  point  to  illustrate  the  degree  of  asymmetry  found  in  the  experiments. 
Compare  the  values  corresponding  to  the  right  edge  of  the  gate  (pressure 
cells  1  through  3)  with  those  corresponding  to  the  left  edge  of  the  gate 
(pressure  cells  7  through  9).  Plate  18  defines  the  curvilinear  coordinate  (s) 
measured  along  the  reference  longitudinal  arcs. 


Table  2 

Averaae  Experimental  Uostream  Pressure  Field  on  Gate  (100-s  test) 

Point 

X,  in. 

S,  in. 

Pressure,  psf 

1 

1.04 

3.52 

21.78 

2 

1.04 

6.40 

32.75 

3 

1.04 

11.60 

56.24 

4 

11.60 

3.52 

20.32 

5 

11.60 

6.40 

38.11 

6 

11.60 

11.60 

58.76 

7 

22.88 

3.52 

24.25 

8 

22.88 

6.40 

36.04 

9 

22.88 

11.60 

62.26 

Note:  Fullv  closed  gat 

The  required  sampling  rate  and  record  length  were  determined  as  follows 
(Bendat  and  Piersol  1993,  Newland  1984).  From  preliminary  eigenproblem 
solution  for  the  gate  structure,  it  was  estimated  that  the  ft-equency  range  of 
interest  goes  as  far  as  70  Hz  to  include  the  first  five  modes  of  vibration. 
Therefore,  the  required  sampling  interval  would  be  1/280  s  for  a  Nyquist  fre¬ 
quency  of  140  Hz.  An  acceptable  accuracy  was  established  by  defining  the 
coefficient  of  variation  of  the  parameter  being  measured  as  C.O.V.  =  u/m  = 
0.14,  which  corresponds  to  the  coefficient  of  variation  of  the  random  vari¬ 
able  Xk  with  100  degrees  of  freedom  (/c  =  100).  Now,  by  estimating  the 
frequency  resolution  as  1  Hz,  the  effective  bandwidth  of  the  spectral  window 
is  given  by  BJ  =  51.  And,  by  defining  B^  =  1/2  Hz,  the  required  record 
length  is  T  =  102  s.  It  was  decided  to  sample  at  a  rate  of  500  sps  for  T  = 
100  s.  Later,  during  the  development  of  the  experimental  project,  even  longer 
records  were  obtained;  T  =  200  s. 
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During  the  first  phase  of  experiments  with  the  physical  model,  the  maxi¬ 
mum  reaction  response  was  obtained  when  a  unit  discharge  of  water  of  50  cfs 
runs  over  the  system  of  fully  raised  gates  with  the  tail  water  pool  maintained  at 
95  ft  (prototype  level)  and  the  headwater  pool  growing  monotonically  from 
about  114  ft  to  119  ft  (prototype  level).  Plate  19  shows  the  minimum  and 
maximum  values  of  the  experimental  pressure  acting  on  the  upstream  side  of 
the  gate  during  the  controlling  100-s  test  (T1080).  The  pressure  field  is 
described  along  three  different  gate  arcs:  the  right  edge  (described  by  a 
sequence  of  dots),  the  middle  arc  (described  by  a  solid  line),  and  the  left  edge 
(described  by  a  sequence  of  small  circles).  Notice  that  the  curves  represent 
parabolic  polynomial  fits  on  three  data  points,  which  explains  the  intersection 
between  the  upper  and  the  lower  bound  curves  for  the  middle  arc  toward  the 
gate  base.  These  minimum  and  maximum  curves  delimit  narrow  bands,  sug¬ 
gesting  very  little  fluctuation  of  the  pressure  at  a  given  location  during  the 
test.  Notice  also  that  the  theoretical  estimation  of  the  pressure  field  is  an 
excellent  upper  bound  over  most  of  the  domain,  except  on  a  small  region 
toward  the  gate  base.  This  observation  together  with  the  analyses  of  tests 
T1108  and  T1109  below,  suggest  the  proposed  design  pressure  line  AB  in 
Plate  19,  as  per  phase  1  of  the  project  (USAEWES  1994a). 

Next,  this  pressure  line  AB  is  tested  against  the  maximum  net  total  pres¬ 
sure  on  the  upstream  side  of  the  gate,  obtained  as  the  algebraic  summation  of 
the  maximum  and  minimum  pressure  fields  on  the  upstream  and  downstream 
sides  of  the  gate,  respectively.  Plate  20  shows  this  maximum  net  pressure  for 
each  reference  longitude  in  the  gate  (i.e.,  right  edge,  middle  arc,  and  left 
edge),  and  for  the  two  200-s  controlling  tests  T1108  (integral  nappe)  and 
T1109  (ventilated  or  broken  nappe).  In  general,  the  test  with  the  integral 
nappe  controls  for  the  upper  half  of  the  gate,  whereas  the  ventilated  nappe  test 
controls  for  the  deeper  half  of  the  gate.  Notice  that  line  AB  represents  a  rea¬ 
sonable  envelope  for  both  the  theoretical  estimations  and  the  maximum  read¬ 
ings  during  the  experiments,  except  for  some  slightly  unconservative  points 
toward  the  gate  ridge.  The  design  total  pressure  field  under  normal  operating 
conditions  is  then  recommended  as 

/?  =  40  +  3.9  5  [psf]  (28) 

where  s  is  the  curvilinear  coordinate  [in.]  along  the  1/15-scale  gate  arc 
(USAEWES  1994a). 

This  design  pressure  field  was  used  as  input  for  separate  finite  element 
analyses  of  the  1/15-scale  physical  model  and  the  prototype.  The  structural 
geometries  were  construct^  analytically  using  the  preprocessor  PATRAN-3 
(PDA  Engineering  1993)  (Plates  21  and  22  display  isometric  views  of  the 
upstream  and  downstream  sides  of  the  finite  element  grid,  respectively);  the 
finite  element  processes  were  executed  by  ABAQUS  (Hibbitt,  Karlsson  and 
Sorensen  1993,  Zienkiewicz  1977)  on  the  structure  standing  on  elastic  sup¬ 
ports  (Teflon  bearings);  and  the  postprocessing  of  output  information  was 
again  effected  by  PATRAN-3. 
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Plate  23  shows  the  displaced  shape  of  the  1/15-scale  physical  model  sub¬ 
jected  to  the  design  pressure  field,  and  Table  3  reports  the  reaction  forces  act¬ 
ing  on  the  three  bearings.  Notice  that  only  the  results  for  one  edge  support 
are  reported  in  the  table  since  the  results  for  the  other  edge  support  are  identi¬ 
cal,  due  to  symmetry.  The  total  resultant  bearing  reaction  turns  out  with 
50  percent  higher  magnitude  than  experimentally  measured,  as  expected, 
because  of  the  conservatism  in  the  design  envelope  for  pressure.  The  corre¬ 
sponding  force  in  the  actuator  was  f^yj  =  131.4  lb  (USAEWES  1994a). 


Table  3 

Montgomery  Point  Gat 
of  Finite  Element  Anah 

e  1/1 5-Scale  Model  (3  s 
i^sis  Bearing  Reactions,  I 

upports);  Summary 
b 

Support  type: 

Edqe 

Center 

Load  cel): 

#1 

#2 

#1 

#2 

Vertical  component: 

25.32 

-83.40 

46.90 

-91.73 

Horizontal  component: 

41.00 

41.00 

41.90 

41.90 

Load  cell  resultant: 

48.19 

92.93 

62.89 

100.85 

Resultant  bearina  reaction: 

Vertical  component: 

-160.99 

Horizontal  component: 

247.80 

Total  reaction  maanitude: 

295.50 

jlNote:  A  minus  sipn  indicates  compressive  force  on  the  support. 

Plate  24  shows  a  plan  view  of  the  prototype  gate  standing  on  five  sup¬ 
ports.  The  design  reactions  at  these  supports  were  obtained  through  finite  ele¬ 
ment  analysis  of  the  structural  idealization  for  the  prototype  gate  (made  out  of 
structural  steel)  subjected  to  the  design  distribution  of  the  operational  hydro- 
dynamic  pressure  field  properly  scaled  (Xp  =  15)  and  acting  on  the  fully 
raised  gate  (70-deg  elevation). 

The  fully  raised  gate  condition  was  recognized  to  control  the  design  of  the 
bearing  supports  among  all  experiments  on  the  physical  model.  This  point  is 
illustrated  in  Plates  25  through  28.  Plate  25  describes  the  variation  of  the  hor¬ 
izontal  (HE)  and  vertical  (VE)  components  of  the  resultant  reaction  at  sup¬ 
ports  1  and  2  with  the  gate  elevation  (X),  and  the  corresponding  engineering 
envelopes  for  design  (notice  that  the  plate  was  actually  recorded  for  bearing 
block  2  in  the  model,  for  the  controlling  test  T1027).  Apparently,  the  maxi¬ 
mum  response  is  obtained  at  the  end  of  the  experiment  when  the  gate  is  fully 
raised.  Plates  26  and  27  lead  to  the  same  conclusion  using  similar  diagrams 
for  the  reaction  at  the  central  block  and  for  the  resultant  reaction  at  supports  4 
and  5  in  the  prototype  gate,  respectively.  Finally,  Plate  28  confirms  the  con¬ 
trolling  fully  raised  gate  position  in  terms  of  the  total  hydraulic  moment  on  the 
torque  tube. 
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Plate  29  displays  an  isometric  view  of  the  upstream  side  of  the  finite  ele¬ 
ment  grid  for  the  prototype,  standing  on  rigid  supports,  and  Table  4  reports 
the  reactive  forces  required  for  the  actual  design  of  the  bearings.  These  reac¬ 
tions  are  for  ideal  conditions  and  do  not  consider  any  misalignment  in  the 
torque-tubes  or  imperfections  in  the  construction.  Correspondingly,  the  force 
in  the  actuator  cylinder  was  F^yj  =  1,699.4  kips  (USAEWES  1994a). 

Plates  30  and  31  show  iso-principal-stress  regions  on  the  gate.  Plate  30 
exhibits  tensile-stress  regions  on  the  upstream  face  of  the  gate,  while  Plate  31 
exhibits  compressive-stress  regions  on  the  elements  of  the  downstream  face  of 
the  gate.  No  indication  of  local  yielding  is  detected  from  the  analysis. 


Table  4 

Montgomery  Point  Gate  Prototype  (5  supports):  Summary  of  Finite 
Element  Analysis  Bearing  Reactions,  kips 

Support  tvpe: 

External  Edae 

Internal  Edqe 

Center 

Vertical  component: 

-200.87 

-221.61 

-376.19 

Horizontal  component: 

388.00 

444.29 

1.259.54 

Resultant  support  reaction: 

436.91 

496.49 

1.314.52 

I  Resultant  bearina  reaction: 

Vertical  component: 

-1,221.15 

Horizontal  component: 

2,924.12 

Total  reaction  maanitude: 

3.168.86 

|lNote:  A  minus  sian  indicates  compressive  force  on  the  support. 
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Plate  14.  Montgomery  Point  Rigid  Gate  Model.  Typical  load  cell  and  corresponding  bridge  wiring 
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Plate  1 5.  Montgomery  Point  Rigid  Gate  Model.  Accelerometer  locations 
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Plate  16.  Montgomery  Point  Rigid  Gate  Model.  Pressure  cell  locations  on  the  upstream  side  of  gate 
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Plate  18.  Montgomery  Point  Rigid  Gate  Model.  Curvilinear  coordinate  (s)  along  the  reference  longitudinal  arcs  (right  edge,  middle 
arc,  and  left  edge)  on  the  upstream  side  of  gate 


EXPERIMENTAL  U/S  PRESSURE  FIELD 
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Plate  19.  Montgomery  Point  Rigid  Gate  Model.  Comparison  of  recommended  design  pressure  field  with  the  theoretical  prediction  and 
its  1 00-s  test  experimental  counterpart 


MAXIMUM  NET  PRESSURE  FIELD  ON  GATE 
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Plate  20.  Montgomery  Point  Rigid  Gate  Model.  Comparison  of  recommended  design  pressure  field  with  the  maximum  net  pressure 
field  on  the  gate  in  200-s  experiments  T1 108  (integral  nappe)  and  T1 109  (ventilated  nappe) 


Scale  Model  F.E.  Grid: 


Plate  21 .  Montgomery  Point  Rigid  Gate  1  /1 5-Scale  Model.  Isometric  view  of  the  upstream  side  of  the  finite  element  grid 
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subjected  to  the  recommended  design  pressure  field 
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Plate  24.  Montgomery  Point  Prototype  Gate.  Plan  view  of  gate  and  location  of  recommended  five  supports 


FORCE,  kips 


TYPE  1 

MONTGOMERY  POINT.  TEST  1027.  BEARING  BLOCK  2 
UNIT  DISCHARGE  -  50.0  cfs 

TAILWATER  ELEVATION  -  95.0  ft 

GATE  OPERATION 

LEFT  GATE  UP 
RIGHT  GATE  UP 

CENTER  GATE  ROTATIING  UP  70  DEGREES  IN  8  MIN. 
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Plate  25.  Montgomery  Point  Prototype  Gate.  Envelopes  for  horizontal  (HE)  and  vertical 
(VE)  components  of  the  resultant  reaction  at  supports  1  and  2  as  a  function  of 
the  gate  elevation  (X).  The  envelopes  are  derived  from  the  records  for  bearing 
block  2  in  the  physical  model  during  the  controling  test  T1027 
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Plate  26.  Montgomery  Point  Prototype  Gate.  Envelopes  for  horizontal  (HC)  and  vertical 

(VC)  components  of  the  resultant  reaction  at  the  central  support  3  as  a  function 
of  the  gate  elevation  (X).  The  envelopes  are  derived  from  the  records  for  bear¬ 
ing  block  4  in  the  physical  model  during  the  controling  test  T1027 
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U/S  POOL  ELEVATION.  FT.  WATER 


FORCE,  kips 


TYPE  1 

MONTGOMERY  POINT.  TEST  1027,  BEARING  BLOCK  6 
UNIT  DISCHARGE  -  50.0  cfs 

TAILWATER  ELEVATION  -  95.0  ft 

GATE  OPERATION 

LEFT  GATE  UP 
RIGHT  GATE  UP 

CENTER  GATE  ROTATIING  UP  70  DEGREES  IN  B  MIN. 


14.0 


28.0 

GATE  ANGLE-E 


42.0 

DEGREES 


56.0 


70.0 


Plate  27.  Montgomery  Point  Prototype  Gate.  Envelopes  for  horizontal  (HE)  and  vertical 
(VE)  components  of  the  resultant  reaction  at  supports  4  and  5  as  a  function  of 
the  gate  elevation  (X).  The  envelopes  are  derived  from  the  records  for  bearing 
block  6  in  the  physical  model  during  the  controling  test  T1027 
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TYPE  1 

MONTGOMERY  POINT.  TEST  1027 
UNIT  DISCHARGE  -  50.0  cfs 

TAILWATER  ELEVATION  -  95.0  ft 

gate  OPERATION 

LEFT  GATE  UP 

cI^^GATE^ROTATIING  up  70  DEGREES  IN  8  MIN. 


3  -  Pool  Elevation 
6  -  Hydraulic  Moment 


M(X)  .  0.0295.x3  _  3  5593.x2  ^  151.0204I.X  -  250, 


“1250.0- 


2^750.0 


28. 0 

gate  ANGLE-2 


42.0 

DEGREES 


Plate  28.  Montgomery  Point  Prototype  Gate.  Envelope  for  the  total  hydraulic  moment  |M) 
as  a  function  of  the  gate  elevation  (X).  The  envelope  is  derived  from  the 
records  for  the  hydraulic  moment  in  the  physical  model  during  the  controling 

test  T1 027 


60 


Chapter  3  Experimental  and  Finite  Element  Analyses 


U/S  POOL  ELEVATION.  FT.  HATER 


16889  Nodes 
5884  Elements 


Plate  29.  Montgomery  Point  Prototype  Gate.  Isometric  view  of  upstream  side  of  finite  element  grid 
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Plate  30.  Montgomery  Point  Prototype  Gate.  Iso-principal-stress  (tensile)  regions  on  the  upstream  face  of  the  gate 
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Plate  31 .  Montgomery  Point  Prototype  Gate.  Iso-principal-stress  (compressive)  regions  on  the  elements  of  the  downstream  face  of  the 
gate 


4  Extreme-Value  Analysis  of 
Experimental  Pressure  Data 


The  objective  of  an  extreme-value  analysis  of  the  experimental  pressure 
data  is  to  provide  the  designer  with  a  rational  estimation  of  the  largest- 
magnitude  hydraulic  loads  to  act  on  the  gate  during  its  useful  life.  Three 
models  are  considered  in  this  study.  The  first  two  models  are  probabilistic  in 
nature  and  are  based  on  well  established  extrapolation  techniques  mostly 
applied  in  offshore  (Cartwright  and  Longuet-Higgins  1956,  Clough  and  Pen- 
zien  1975)  and  wind  engineering  (Davenport  1963,  Clough  and  Penzien  1975, 
Simiu  and  Scanlan  1986),  respectively.  The  third  model  aims  at  the  determi¬ 
nation  of  deterministic  design  pressure  fields  involving  the  results  of  many 
tests  subjecting  the  gate  physical  model  to  multiple  unusual  and  extreme  load¬ 
ing  conditions  during  the  second  phase  of  the  experimental  project. 


Probability  Distribution  for  Maxima 


Assuming  a  zero-mean  stationary  Gaussian  process  $(t)  (Benjamin  and  Cor¬ 
nell  1970,  Crandall  and  Mark  1963,  Naudascher  1991)  having  an  arbitrary 
spectral  density  Sj(f),  Cartwright  and  Longuet-Higgins  (1956)  have  derived 
the  probability  density  function  for  maxima  in  the  process  as 


p(v)  = 


1 

(2u)‘/2 


7)(l-€^)''^ 


-00 


(29a) 


where  the  nondimensional  random  field  is  defined  by 


V  = 


1/2 


(29b) 
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the  parameter  e  is  given  by 


/Ho  7^4 


(29c) 


and  /H,-,  i  =  0,1... 4,  are  the  moments  of  the  spectral  density  function  about 
the  origin  in  the  frequency  domain. 

Plate  32  shows  this  probability  density  function  for  different  values  of  the 
parameter  e.  It  can  be  shown  that  the  parameter  e  lies  in  the  range  (0,1). 

For  the  limiting  case  e  -*0,  this  density  function  reduces  to  the  form  of  a  Ray¬ 
leigh  distribution,  whereas  for  the  limiting  case  e  ^1,  the  resulting  distribution 
of  maxima  approaches  the  form  of  a  Gaussian  distribution  (Cartwright  and 
Longuet-Higgins  1956,  Clough  and  Penzien  1975). 

Test  T1080  was  selected  as  a  pilot  experiment  to  determine  the  nature  of 
the  parameter  e  in  the  context  of  the  experimental  project.  The  process  is 
assumed  stationary  and  ergodic  (Yang  1986,  Crandall  and  Mark  1963).  The 
time-average  (mean)  for  each  pressure-cell  record  in  the  test  was  first  com¬ 
puted  and  extracted  from  the  data  (refer  to  the  third  column  in  Table  5).  The 
power  spectral  density  function  was  then  generated  for  the  upstream,  down¬ 
stream,  and  net  upstream  pressure  records,  respectively,  for  each  cell  on  the 
instrumented  gate.  Plates  33  to  35  show  representative  spectral  densities— in 
this  particular  case,  for  pressure  cell  4.  The  area  under  each  spectral  density 
(ttiq),  as  well  as  the  second  and  fourth  moments  (m2  and  m^),  were  calculated 
by  means  of  a  computer  program  under  MATLAB  (MATHWORKS  1993) 
(Table  5). 

The  last  column  of  Table  5  shows  the  resulting  value  of  the  parameter  e  for 
each  pressure  cell  recorded  during  the  test.  Notice  that  most  results  are  above 
0.9.  Thus,  for  all  practical  purposes,  the  energy-content  parameter  e  may  be 
assumed  as  one  (e  =  1)  for  the  random  field  under  study  (17),  and  the  distribu¬ 
tion  for  the  peaks  are  assumed  to  be  Gaussian.  Therefore,  the  maxima  for  the 
pressure  field  may  be  predicted  using  standard  tables  for  the  normal  distribu¬ 
tion,  for  any  confidence  level  selected  by  the  designer.  For  example,  based 
on  test  T1080,  the  maximum  value  predicted  for  the  net  upstream  pressure  on 
cell  1  is  its  mean  value  plus  1  -Vm^,  2  -Vm^,  or  3  -V/hq  ,  with  a  confi¬ 
dence  level  of  68.27,  95.45,  or  99.73  percent,  respectively  (Benjamin  and 
Cornell  1970,  Simiu  and  Scanlan  1986).  Assuming  a  95.45-percent  confi- 
dence  level,  this  maximum  net  upstream  pressure  is  8.3368  +  2  •  V4.8196 
in.  of  water  =  12.73  in.  of  water,  or,  approximately,  66.2  psf.  The  exact 
95.45-percent  confidence-level  extreme-pressure  would  be  given  by  8.3368  -t- 
17V4.8I96  in.  of  water,  where  1/*  is  the  abscissa  in  Plate  32  for  the  curve  cor¬ 
responding  to  €  =  0.9122  such  that  the  area  under  the  curve  from  -00  to  17*  is 
97.75  percent.  This  technique  was  verified  to  give  results  close  to  those 
provided  by  the  second  model,  and  since  the  second  model  (described  next)  is 
computationally  simpler,  the  present  model  was  not  elaborated  any  further. 
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Table  5 

Computation  of  Energy  Parameter  e  for 
Duration:  100  s) 


Test  T1 080  (Sampling  Rate:  500  sps; 
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Probability  Distribution  for  Extreme  Values 


The  second  model  is  the  preferred  probabilistic  extrapolation  method  for 
being  computationally  efficient  and,  to  a  certain  extent,  experimentally  veri¬ 
fied  in  the  context  of  wind  engineering.  The  model  was  originally  derived  by 
Davenport  (1963)  and  consists  of  a  modification  of  the  Gumbel  Type  I 
extreme-value  distribution  of  the  largest  values  (Benjamin  and  Cornell  1970, 
Clough  and  Penzien  1975).  In  this  method,  it  is  assumed  that  the  parameters 
of  the  distribution  can  be  obtained  by  replacing  the  mean  and  the  mean  square 
value  of  the  random  variable  pressure  (X)  by  the  corresponding  statistics  of 
the  sample.  Specifically, 


a  = 


(30a) 


and 


{i  =  X  -  0.5772  a 


(30b) 


where  X  and  s  are  the  sample  mean  and  the  sample  standard  deviation, 
respectively. 

With  these  estimators,  the  following  estimator  can  be  constructed 


Gj^ip)  =  50  -  0.5772)  ^ 

TT 


(31a) 


where 

y=  -ln(-ln/j)  ,  (31b) 

and  p  is  the  probability  of  exceeding  the  target  pressure  level  within  a  prede¬ 
fined  subinterval  of  the  mean  recurrence  interval. 

Under  the  assumption  that  the  random  variables  X  and  s  are  asymptotically 
normally  distributed,  it  can  be  shown  that,  for  large  samples  of  size  n,  the 
standard  deviation  of  the  estimator  in  Equation  31  is  given  by 


SD 


—  +  1.13960  -  0.5772)—  + 

6  ye 


1.10  -  0.5772)2 


1/2 


a 


(32) 
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Then  the  approximate  statement  can  be  made  that  the  interval  0^^)  +  SD 
corresponds  to  the  68-percent  confidence  level,  ±  2'SD  corresponds  to 
the  95-percent  confidence  level,  and  ±  3  SD  corresponds  to  the 
99.7-percent  confidence  level. 

For  example,  considering  the  controlling  test  T1 108  (sampling  rate: 

500  sps;  duration:  200  s),  it  was  decided  to  envision  the  records  for  pressure 
as  composed  by  series  of  2-s  lapses.  So,  the  starting  objective  is  to  generate  a 
table  of  two-secondly  maximum  total  hydraulic  pressures.  Since  each  record 
contains  100,000  data  points,  and  each  2-s  lapse  contains  1,000  points,  then 
the  table  of  two-secondly  maxima  contains  1()0  data  points  from  which  the 
sample  statistics  (mean  X  and  standard  deviation  s)  are  computed. 

If  it  is  decided  to  design  the  gate  such  that  the  allowable  stress  under  total 
hydraulic  pressure  is  exceeded  with  a  50-year  mean  recurrence  period,  i.e., 

50  X  365  X  24  X  1,800  =  7.884  £-1-8  two-s  recurrence  interval,  then  the 
probability  of  exceedance  in  any  2-s  subinterval  is 


P  = 


1 

7.884  £h8 


1.2684  B-9 


The  probability  that  the  allowable  stress  be  attained  at  least  once  in  a  given 
number  of  years  is  given  by  1  -  (1  -p)  (1*365x24x1, soo)  appears  listed  in 

Table  6  (Benjamin  and  Cornell  1970). 


Table  6 

Probability  that  Allowable  Stress  be  Exceeded  at  Least  Once  in  a 
Given  Time  Interval,  50-year  Mean  Recurrence  Interval 

Time  Interval,  Years 

Probability  of  Exceedance,  Percent 

1 

2.0 

5 

9.5 

10 

18.02 

20 

32.79 

25 

39.15 

50 

63.0 

100 

86.3 

A  computer  program  under  MATLAB  (MATHWORKS  1993)  was  devel¬ 
oped  for  the  automatic  execution  of  Equations  30  through  32.  Plate  36  shows 
the  underlying  probability  distribution  for  the  extreme  value  of  the  net  pres¬ 
sure  field  on  the  upstream  face  of  the  physical  model  at  the  location  of  cell  1, 
as  an  example.  Plates  37  and  38  show  the  prediction  of  this  extrapolating 
mathematical  model  for  the  same  pressure  field  on  the  upstream  side  of  the 
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gate  for  the  pilot  test  T1080  and  for  the  controlling  test  T1108,  respectively. 
Each  plate  presents  the  extreme-pressure  curves  for  the  reference  longitudinal 
arcs  of  the  gate— i.e.,  the  right  edge  (sequences  of  dots),  the  middle  arc  (solid 
lines),  and  the  left  edge  (sequences  of  circles),  looking  in  the  downstream 
direction  into  the  upstream  side  of  the  gate  (Plate  18).  The  curves  denoted  as 
1,  2,  and  3  represent  parabolic  polynomial  fits  on  three  data  points  for  the  68, 
95,  and  99.7-percent  confidence  levels,  respectively. 

The  plots  in  Plates  37  and  38  definitely  show  similar  trends. 


Pressure  Envelopes— Phase  2  of  the  Experimental 
Project 

A  second  phase  of  the  experimental  project  was  devised  to  further  examine 
the  response  of  the  gate  under  extreme  conditions  of  discharge  and  pool  eleva¬ 
tions,  and  unusual  operation,  or  both  (USAEWES  1994b).  In  particular,  dur¬ 
ing  the  experiments  in  the  second  phase,  the  unit  discharge  was  generally 
variable  and  larger  than  50  cfs  (the  unit  discharge  was  maintained  fixed  during 
the  first  phase).  Also,  in  this  second  phase  of  the  project,  conditions  involv¬ 
ing  extensive  motion  over  the  whole  range  of  inclinations  of  the  instrumented 
gate  were  considered. 

For  each  of  37  experiments  in  the  second  phase,  a  plot  presenting  the  most 
relevant  deterministic  information  on  the  pressure  field  was  generated 
(USAEWES  1994b).  However,  in  reviewing  the  data  from  the  Phase  2  study, 
the  WES  engineering  team  discovered  that  a  key  pressure  cell  had  malfunc¬ 
tioned  during  the  experiments.  The  controlling  tests  of  the  Phase  2  study 
(Table  7)  were  repeated  and  new  design  pressure  fields  on  the  gate  were 
generated. 

The  information  summarized  next  complements  the  pressure  contours  pre¬ 
sented  in  Report  1  of  this  study  (Fletcher  and  de  Bejar  1995).  As  before, 
each  plot  represents  the  most  relevant  deterministic  information  for  the  indi¬ 
vidual  tests  in  the  series.  In  general,  for  each  reference  longitudinal  arc  in  the 
gate  (i.e.,  the  right  edge,  the  middle  arc,  and  the  left  edge),  a  parabolic 
polynomial  has  been  fitted  to  the  three-point  spectral  data  extracted  from  the 
net  pressure-cell  time-histories  of  response  (duration:  100  s  per  test).  Plates 
39  and  40  present  the  results  for  the  controlling  individual  tests  in  the  series* 
(analytical  expressions  for  each  design  pressure  line  appear  at  the  bottom  of 
each  plate).  Notice  that,  although  pressure  field  results  are  now  smaller  in 
magnitude  than  previously  reported  (USAEWES  1994b, c),  the  tests  controlling 
the  design  remain  the  same,  i.e.,  test  T1 192  (before  test  T1 160)  for  the  upper 
region  of  the  gate  and  tests  T1 179  and  T1 180  (before  tests  T1 139  and  T1 140) 
for  the  intermediate  and  deep  regions  of  the  gate,  respectively. 


Report  1  of  this  study  describes  the  results  for  all  tests  in  Phase  2  of  the  study. 
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Table  7 

Montgomery  Point  Gate  Model,  Phase  2,  New  Tests  to  Replace 

Some  Defective  Pressure-Cell  Records 

Upper 

Lower 

Initial 

Pool 

Pool 

Center 

Center 

Test 

Previous 

q’ 

Eiev., 

Elev. 

Left  Gate 

Gate 

Gate 

Right  Gate 

No. 

Test  No. 

cfs/ft 

ft 

ft 

Position 

Position 

Action 

Position 

1179 

1139 

53.3 

119 

95 

Raised 

Raised 

Hold 

Raised 

1180 

1140 

111.1 

119 

95 

Raised 

Raised 

Hold 

Lowered 

1181 

1141 

144 

116.8 

95 

Lowered 

Raised 

Hold 

Lowered 

1183 

1143 

Varies 

119 

95 

Raised 

Raised 

Lowering 

Raised 

1184 

1144 

Varies 

117 

95 

Raised 

Raised 

Lowering 

Lowered 

1187 

1145 

Varies 

116.8 

95 

Lowered 

Raised 

Lowering 

Lowered 

1188 

1152 

Varies 

115 

95 

Raised 

Lowered 

Raising 

Raised 

1189 

1153 

Varies 

115 

95 

Raised 

Lowered 

Raising 

Lowered 

1190 

1154 

Varies 

119 

95 

Raised 

Lowered 

Raising 

Raised 

1191 

1155 

Varies 

117 

95 

Raised 

Lowered 

Raising 

Lowered 

1192 

1160 

Varies 

121 

115 

Raised 

Raised 

Lowering 

Raised 

^  Unit  discharge  based  on  a  90-ft-wide  approach  channel. 

All  these  tests  feature  the  central  and  left  gates  in  raised  position.  In  tests 
1179  and  1192,  the  right  gate  is  also  held  raised;  but  in  the  case  of  test  1179, 
the  central  gate  (instrumented)  is  kept  stationary,  whereas  in  test  1 192,  the 
central  gate  goes  to  a  lowered  position  during  the  test  (with  the  upstream  pool 
maintained  above  the  normal  level,  simulating  improper  operation).  In  test 
1180,  the  right  gate  is  held  steadily  in  lowered  position  (with  the  upstream 
pool  maintained  at  the  peak  level,  simulating  extreme  normal  operation). 


The  envelopes  in  Plates  40  and  41  refer  to  the  prototype  and  are  provided 
as  information  to  the  designer.  They  are  expressed  algebraically  as  follows: 

For  the  Right  Edge: 

-80.0  ^  +  1,000  [psf]  0  <  5  <  5  ft 


-1-18.75  ^ -I-  506.25  [psf|  5  ft  <  5  <  13  ft 

-1-9.649  P-  -  169.1228  ^  -I-  1,317.895  [psf]  13  ft  <  ^  <  20  ft 
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For  the  Middle  Arc: 


+  16.2  P-  -  348.5  J  +  1,996.2  [psf]  0  <  ^  <  4.5  ft 

+  11.533  5  +  704.1  [psf]  4.5  ft  <  5  <  15  ft 

+  10.9  -  193  5  +  1,319.6  [psf]  15  ft  <  5  <  17.5  ft 

For  the  Left  Edge: 

-100.0  J  +  1,500  [psf]  0  ^  J  <  7  ft 

+  800  [psf]  7  ft  <  ^  <  14.5  ft 

+  127.27  j  -  1,045.46  [psf]  14.5  ft  <  ^  <  20  ft, 

where  s  is  the  depth  of  the  point  in  question  measured  along  the  gate  arc. 

These  input  forcing  functions  should  be  used  in  the  context  of  Allowable 
Stress  Design  for  structural  steel,  i.e.,  there  is  no  need  for  further  factoring 
the  loads  since  sufficient  conservatism  has  already  been  built  into  the  estima¬ 
tion  procedure  (AISC  1986). 

Notice  that  the  pressure  levels  as  determined  in  Phase  2  of  the  project  take 
precedence  over  all  results  from  Phase  1,  even  over  the  results  of  the  extreme- 
value  analysis.  Larger  unit  discharges  are  considered  responsible  for  this  sub¬ 
stantial  increase  in  the  magnitude  of  the  recommended  design  pressure  field  on 
the  gate. 
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Plate  32.  Probability  density  function  for  the  peaks  of  the  nondimensional  random  field  /?  (/?  =  for  different  values  of  the 

parameter  e 


-XP 


Spectral  Density 
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Plate  34.  Spectral  density  for  the  pressure  at  cell  4  on  the  downstream  face  of  the  gate  model  during  test  T1080 
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Plate  35.  Spectral  density  for  the  net  pressure  on  the  upstream  face  of  the  gate  model  at  the  location  of  cell  4  during  test  T1080 


Plate  36.  Probabilistic  description  of  the  extreme  value  of  the  net  pressure  on  the 

upstream  face  of  the  gate  model  at  the  location  of  cell  1  during  test  T1 108  [in. 
of  H2O].  la)  Probability  density  function  IGumbel  Type  I,  of  the  largest  values), 
and  (b)  the  corresponding  probability  distribution  function 
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Plate  37.  Extreme  pressure  field  as  predicted  by  test  T1080  on  the  three  reference  longi¬ 
tudinal  arcs  on  the  gate  skin,  i.e.,  the  right  edge,  the  middle  arc,  and  the  left 
edge.  Level  1 :  68-percent  confidence  bound;  level  2:  96-percent  confidence 
bound;  and  level  3:  99.7-percent  confidence  bound 
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Total  PRESSURE  [psf] 


Plate  38.  Extreme  pressure  field  as  predicted  by  test  T1 108  on  the  three  reference  longi¬ 
tudinal  arcs  on  the  gate  skin,  i.e.,  the  right  edge,  the  middle  arc,  and  the  left 
edge.  Level  1;  68-percent  confidence  bound;  level  2:  96-percent  confidence 
bound;  and  level  3;  99.7-percent  confidence  bound 
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DESIGN  PRESSURE  LINE  AS  PER  TEST  T1 179  (former  T1 1 39)  DESIGN  PRESSURE  LINE  AS  PER  TEST  T1 180  (former  T1 140) 
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Plate  39.  Montgomery  Point  Prototype  Gate.  Design  pressure  line  as  per  tests  T1 1 79  and  T1 1 80 


DESIGN  PRESSURE  LINE  AS  PER  TEST  T1192  (former  T1 160)  DESIGN  PRESSURE  LINES  FOR  MIDDLE  ARC 
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Plate  40.  Montgomery  Point  Prototype  Gate.  Design  pressure  line  as  per  test  Til 92.  Pressure  envelope  for  middle  arc 


DESIGN  PRESSURE  LINES  FOR  LEFT  EDGE  DESIGN  PRESSURE  LINES  FOR  RIGHT  EDGE 
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Plate  41 .  Montgomery  Point  Prototype  Gate.  Pressure  envelope  for  left  edge.  Pressure  envelope  for  right  edge 


5  Gate  Eigenproperties 


In  this  chapter,  the  experimental  identification  of  natural  frequencies  of 
vibration  and  corresponding  mode  shapes  is  described.  Knowledge  of  these 
structural  properties  for  both  the  gate  model  and  the  prototype  system  is 
important  to  verify  parameters  entering  future  modal  analyses  of  the  gate;  to 
compare  these  parameters  in  relation  to  the  actual  prototype  gate  as  it  is  actu¬ 
ally  constructed;  and  to  keep  track  of  structural  deterioration  over  the  life  of 
the  gate  in  service,  due  to  weathering,  corrosion,  local  fatigue  of  component 
elements,  or  simply  due  to  normal  wear  (safety  assessment)  (Crandall  and 
Mark  1963). 

A  variety  of  mutually  checking  methods  were  implemented  in  this  exer¬ 
cise.  A  finite  element  grid  for  the  physical  model  standing  on  Teflon  bearings 
on  load-cell  foundations  was  first  used  to  identify  its  first  10  modes  of  vibra¬ 
tion  (Zienkiewicz  1977;  Clough  and  Penzien  1975;  Hibbitt,  Karlsson  and 
Sorensen  1993).  Plates  42  through  51  describe  the  results  for  each  relevant 
mode  in  sequence.  Plates  52  through  59  similarly  describe  the  eigensolution 
for  the  prototype  gate  (structural  steel)  standing  on  a  rigid  foundation,  for  the 
first  seven  modes  in  sequence.  As  expected,  the  model  exhibits  much  more 
rigid  characteristics  than  the  prototype,  but  the  same  general  mode-shapes. 

A  dry  Twang  test  (pull-back  test)  (Clough  and  Penzien  1975)  was  also 
conducted  on  the  physical  model  (T0137).  Plate  60  shows  the  time-history  of 
a  typical  parameter  of  response,  say,  the  bending  moment  at  the  root  of  the 
right  operating  arm  (looking  in  the  downstream  direction),  in  its  free-vibration 
decay.  Plate  61  shows  the  spectral  density  for  the  signal.  The  location  of  the 
first  peak  of  this  function  identifies  the  fundamental  frequency  of  vibration  of 
the  model  (~  46.875  Hz). 

The  natural  properties  of  the  physical  model  were  also  identified  using  a 
scanning  Laser  vibrometry  system  (LAZON,  by  Zonic  Co.).  Plate  62a  shows 
the  test  setup  (ZONIC  +  AND  1993).  The  instrumented  (central)  gate  is 
excited  by  a  variable-frequency  electro-mechanical  shaker  connected  to  its 
downstream  face,  while  its  upstream  skin  is  scanned  by  the  Laser  beam  of  a 
vibrometer  during  its  dynamic  response  (Plate  62b).  Plate  63,  for  example, 
shows  the  frequency-response  function  for  velocity  of  several  previously 
defined  points  on  the  gate  skin  (nodes  1,  7,  10,  and  13  in  Plate  64),  whereas 
Plate  65  shows  an  expansion  of  the  function  in  the  important  low-frequency 
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part  of  the  domain  (up  to  120  Hz),  for  nodal  points  19  and  21.  The  identified 
eigenvalues  appear  listed  in  Table  8.  The  first  two  natural  frequencies  are  in 
close  agreement  as  predicted  by  the  various  techniques.  Plates  66  and  67 
illustrate  the  corresponding  mode-shapes.  Plate  68  is  an  example  of  animation 
of  a  higher  mode. 


Table  8 

Natural  Frequencies  for  the  Physical  Model,  Hz 

Mode 

Finite  Element 

Twang  Teet 
DRY 

Accelerometer 

LAZON 
0-100  Hz 

Description 

1 

46.91 

46.875 

50.55 

49.66 

Rigid-Plate  Rocking 

(1) 

- 

- 

55.5 

56.99 

Rigid-Plate  Rocking 

2 

74.66 

- 

79.8 

81.07 

Rigid-Plate  Twist 

(2) 

- 

- 

83.96 

85.45 

Rigid-Plate  Twist 

3 

137.96 

- 

103.66 

102.07 

Rigid-Plate  Twist 

- 

- 

- 

175.3 

- 

Bending 

- 

- 

- 

204.2 

- 

Torsion 

Plates  69  and  70  show  very  little  variation  in  the  frequency-response 
function  from  the  situation  in  which  all  gates  are  raised  to  the  condition  in 
which  only  one  lateral  gate  is  down.  However,  Plate  71  shows  strong  energy- 
suppression  in  the  region  of  the  frequency-domain  beyond  about  60  Hz  for  the 
condition  in  which  both  lateral  gates  are  down  and  only  the  instrumented  gate 
is  raised.  Plate  72  shows  the  negligible  difference  in  the  frequency-response 
function  between  the  instrumented  gate  with  and  without  the  trash  rack. 

Plates  73  through  78  illustrate  the  identification  of  the  fundamental 
frequency  of  free  vibrations  for  the  gate  physical  model  in  the  upstream-sub¬ 
merged  condition  from  accelerometer  recordings  in  Twang  test  TO  139  (fj  « 

39  Hz).  The  system  becomes  more  flexible  (compare  with  fj  »  47  Hz  in 
dry)  due  to  the  “added  mass”  effect  of  the  upstream  water. 

The  identification  of  the  equivalent  viscous  damping  ratio  for  both  the  dry 
and  upstream  submerged  conditions  is  discussed  in  Chapter  7  in  the  context  of 
the  engineering  model,  and  the  analysis  of  forced  vibrations  of  the  gate  in 
operation  is  considered  in  Appendix  A  (Spectral  Analysis). 
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}formed_plot:  LC=6.7-RES=1.1-P3/PATRAN  R.1.2-Deforniation-ABAQUS-11 -Aug-94  13:4 
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Plate  42.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  1 
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Plate  43.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  2 


C=6.9-RES 


Plate  44.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  3 


jformedjjlot;  LC«6.10*RES=1.1-P3/PATRAN  R.1.2-Deformation-ABAQUS-11 -Aug-94  09: 
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Plate  45.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  4 


jformed _plot;  LC»6.11-RES=1.1*P3/PATRAN  R.1.2-Deformation-ABAQUS-11-Aug-94  10; 
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Plate  46.  Finite  element  eigensolution  for  1/15*scale  model.  Relevant  mode  of  vibration  5 
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Plate  47.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  6 


jformed _plot:  LC«*6.13-RES=1.1-P3/PATRAN  R.1.2-Deformatlon-ABAQUS- 11 -Aug-94  11: 
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Plate  48.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  7 
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Plate  49.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  8 


}formed_plot:  LC=6.15-RES=1.1-P3/PATRAN  R.1.2-Defonnation-ABAQUS-11 -Aug-94  11: 
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Plate  50.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  9 
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Plate  51.  Finite  element  eigensolution  for  1/1 5-scale  model.  Relevant  mode  of  vibration  10 
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Plate  52.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  1 
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Plate  53.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  2 
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Plate  54.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  3 
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Plate  55.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  3 
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Plate  56.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  4 
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Plate  57.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  5 
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Plate  58.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  6 
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Plate  59.  Finite  element  eigensolution  for  prototype.  Relevant  mode  of  vibration  7 
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Plate  61 .  Spectral  density  function  for  the  time-history  of  response  in  Plate  60 
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(b)  Laser  Position 


Plate  62.  Scanning  Laser  Vibrometry  System,  (a)  Test  setup,  and  (b)  Laser  position 
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Plate  63.  Scanning  Laser  Vibrometry  System.  Frequency-response  function  for  selected  nodes  on  the  gate  skin 


Top  Edge  of  Surface 
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Plate  64.  Scanning  Laser  Vibrometry  System.  Node  definition  on  the  gate  skin 
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Plate  65.  Scanning  Laser  Vibrometry  System.  Frequency-response  function  for  the  velocity  of  nodes  19  and  21  on  the  gate  skin 
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Plate  66.  Scanning  Laser  Vibrometry  System.  Mode-shape  1 
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Plate  67,  Scanning  Laser  Vibrometry  System.  Mode-shape  2 


Plate  68.  Scanning  Laser  Vibrometry  System.  Mode-shape  5 


Frequency  response 
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when  all  three  gates  are  raised 


Frequency  response  function 
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Plate  70.  Scanning  Laser  Vibrometry  System.  Frequency-response  function  for  the  acceleration  of  nodes  3  and  4  on  the  gate  skin  when 
one  lateral  gate  is  down  and  the  other  two  are  raised 


Frequency  response 
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Plate  71 .  Scanning  Laser  Vibrometry  System.  Frequency-response  function  for  the  acceleration  of  nodes  4  and  6  on  the  gate  skin 
when  both  lateral  gates  are  down  and  the  instrumented  (central)  gate  is  raised 


Frequency  response  function 
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Plate  72.  Scanning  Laser  Vibrometry  System.  Comparison  between  the  frequency-response  functions  for  the  acceleration  of  nodes  1 
and  2  on  the  gate  skin  with  and  without  a  leaning  trash  rack 


[g] 


Plate  73.  Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Transverse  time 
history  of  response  for  triaxial  accelerometer  1  (Twang  test  T0139) 
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Plate  74.  Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Spectral  density 
for  transverse  time-history  of  response  for  triaxial  accelerometer  1  (Twang  test 
T0139) 
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Plate  75.  Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Transverse  time- 
history  of  response  for  triaxial  accelerometer  6  (Twang  test  T0139) 
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Frequency 


Plate  76.  Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Spectral  density 
for  transverse  time-history  of  response  for  triaxial  accelerometer  6  (Twang  test 
T0139) 


Chapter  5  Gate  Eigenproperties 


121 


Plate  77. 


Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Time-history  of 
response  for  uniaxial  accelerometer  4  (Twang  test  T0139) 
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Plate  78.  Free  vibration  of  gate  with  upstream  side  submerged  in  water.  Spectral  density 
for  time-history  of  response  for  uniaxial  accelerometer  4  (Twang  test  TO  139) 
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6  Engineering  Model 


In  this  chapter,  an  engineering  model  is  derived  to  estimate  the  response 
of  the  gate  system  subjected  to  the  upstream  hydraulic  pressure  field.  This 
model  is  based  on  the  principles  of  structural  dynamics  applied  to  an  idealiza¬ 
tion  of  the  gate  system  conceived  as  a  rigid-body  torque-tube  gate  standing  on 
simple  supports  and  being  acted  upon  by  the  hydraulic  pressure  field  and  by 
the  visco-elastic  system  of  resisting  forces  provided  by  the  operating  arms  and 
the  horizontal  actuator  (Clough  and  Penzien  1975).  The  model  is  general,  but 
the  specific  numerical  applications  illustrating  the  development  refer  to  the 
1/ 15-scale  physical  model. 

The  first  section  presents  the  deterministic  derivation  of  the  steady-state 
component  of  the  gate  response— in  particular,  the  gate  rotation  and  the 
moment  at  the  root  of  the  operating  arms.  The  second  section  presents  a  ran¬ 
dom  vibration  model  for  the  stochastic  description  of  the  fluctuations  of  struc¬ 
tural  response  about  the  steady-state  component  (Naudascher  1991,  Yang 
1986).  Both  models  compare  satisfactorily  with  the  corresponding  experimen¬ 
tal  verifications. 


Deterministic  Mathematical  Model 

Consider  the  rigid-body  gate  in  closed  position  (angle  2a  =  70  deg)  shown 
in  Plate  79,  standing  on  Teflon  bearings  and  being  excited  by  an  external 
dynamic  pressure  field  p  and  by  its  own  weight  W  (mass  M  =  Wig).  The 
overturning  moment  of  the  external  forces  about  the  center  of  the  torque  tube 
is  resisted  by  the  eccentric  force  exerted  by  the  horizontal  actuator  with  stiff¬ 
ness  K  and  equivalent  viscous  damping  C,  acting  at  the  far  end  of  the  operat¬ 
ing  arms  with  flexural  rigidity  El  and  length  L.  The  analytical  model  in 
Plate  80  is  then  proposed.  The  gate  is  idealized  as  a  rigid  circular  arc  (radius 
a)  spanning  from  an  elevation  6^  =  20  deg  at  0  to  an  elevation  ^  = 

72.75  deg  at  D,  and  its  mass  is  lumped  at  the  centroidal  point.  The  torque 
tube  reduces  to  point  O,  representing  the  rigid  connection  of  the  operating 
arms  to  the  base  of  the  gate. 

The  displacements  are  discretized  in  this  model  by  the  three  coordinates 
shown  in  Plate  81,  i.e,,  rotations  0,  and  $2  at  nodes  O  and  B,  respectively. 
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and  the  linear  displacement  A  of  the  end  B  of  the  operating  arms  (Rubinstein 
1966,  Zienkiewicz  1977).  Only  the  strain  energy  associated  with  bending 
deformations  of  the  operating  arms  is  considered  (Boresi  and  Sidebottom 
1985,  Rubinstein  1966). 

The  stiffness  matrix  of  the  system  may  be  formulated  as  (Przemienieki 
1985,  Rubenstein  1966) 


2 


3 

L 


{K\ 


2  El  1 
L 


2 


2 

L 


A  2- 

Z  L  (£2  2EI  ^ 


(33) 


Elimination  of  coordinates  62  and  A  from  the  corresponding  stiffness  for¬ 
mulation  by  static  condensation  (Clough  and  Penzien  1975,  Przemienieki 
1985)  produces  the  generalized  stiffness  coefficient  for  the  degree-of-ffeedom 
only  as 


1 

. _ \ _ 

3Ef  kL^co&l 


(34) 


A  similar  procedure  applied  to  the  viscous  damping  matrix  renders  the  gen¬ 
eralized  viscous  damping  coefficient  associated  with  the  degree-of-freedom 
only  as 


a  =  CL^co^a  (35) 

Plate  82  shows  the  generalized  single  degree-of-fteedom  system  (<^  a  0j) 
vibrating  under  the  action  of  the  vertical  and  horizontal  components  of  the 
resultant  hydraulic  force  on  the  gate  (7?,  and  R2,  respectively),  the  gate  own 
weight  (W),  and  the  resisting  elastic,  viscous,  and  inertial  moments  (Af^, 
and  Af,-,  respectively)  at  the  gate  base.  Notice  in  the  plate  the  viscous  damper 
Cp  representing  the  drag  force  on  the  gate,  opposing  its  velocity  through 
water.  The  degree  of  freedom  </>  is  defined  as  the  rotation  of  the  gate  with 
respect  to  its  static  weight-deflected  configuration  (Clough  and  Penzien  1975). 
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The  dynamic  equilibrium  of  the  rigid-body  gate  is  formulated  by  applying 
the  virtual  work  principle  on  the  structural  system  subjected  to  the  virtual  dis¬ 
placement  5<l>  (a  small  variation  from  the  instantaneous  position  of  equilibrium 
determined  by  <f>)  (Mura  and  Koya  1992).  Table  9  lists  the  forces  on  the  gate, 
the  corresponding  virtual  displacements,  and  their  virtual  work  contributions 
to  the  energy  equation  E  (6  Work)  =  0,  i.e.. 


4  ^  +  (CL-cos^a  +  ^  + 


/ 


1 


1 

KL^cos^a 


-  +  R,r^  -  Ww 


(f)  =  R^r^  +  R,r^  , 


(36a) 


) 


where 


=  mass  moment  of  inertia  of  the  gate  about  the  torque-tube  axis 

rj  and  r-^  =  lever  arms  of  forces  i?i  and  R^  about  the  same  axis, 
respectively 


Table  9 

Dynamic  Equilibrium,  Formulation  of  the  Principle  of  Virtual 
Displacements 

Action 

Virtual  Displacement 

Virtual  Work 

^0 

II 

S<p 

‘/o  •  0  '  ^0 

C*  .<p 

6<p 

-C  cos^  a  .  0 ' 60 

~  ^D*  •  ^ 

6<p 

-Cp*  .  0 ' 60 

=  /c* .  ^  +  w  •  V 

60 

^  If,  -6  (f)  Wv(>i) 

3^/  KL^cos^a 

(r-|  +  r2  0)  '  60 

Rl  {r-]  +  r2  0)  *  60 

R, 

(T2  ‘  0)  '  ^<P 

^2  ^^2  ”  *^1  0^  *  ^0 

w 

(v  +  w  •  0)  *  <50 

W  (v  +  w  •  0)  •  <50 

Recognizing  that  the  second-order  terms  within  the  parenthesis  in 
Equation  36, 
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+  /Zj/i  -  Ww  , 


contribute  to  its  numerical  value  with  a  negligible  amount,  and  that  the 
fraction 


1 

KL^cos^a 


is  also  negligibly  small,  the  equation  of  motion  may  be  reduced  in  this  case  to 


(36b) 


where 


3Er 

loL 


(36c) 


is  the  square  of  the  natural  circular  frequency  of  vibration,  and 


I  = 


CL^co^a  +  Cp 

4 


(36d) 


is  the  overall  equivalent  viscous  damping  ratio. 

Plate  83  shows  a  relevant  portion  of  the  free-vibration  decay  of  the  gate 
during  the  Twang  test  TO  137  in  dry  air  and  the  application  of  the  logarithmic 
decrement  method  for  the  approximate  determination  of  the  equivalent  viscous 
damping  ratio  (^  =  12.0  percent,  on  the  average).  Likewise,  Plate  84  shows 
a  relevant  portion  of  the  free-vibration  decay  of  the  gate  during  the  Twang  test 
T0139,  with  water  on  its  upstream  side  only,  and  the  application  of  the  loga¬ 
rithmic  decrement  method  for  the  approximate  determination  of  the  equivalent 
viscous  damping  ratio  (^  =  9.5  percent,  on  the  average)  (Clough  and  Penzien 
1975). 

The  computation  of  was  effected  with  the  assistance  of  the  preprocessor 
PATRAN-3  (PDA  Engineering  1993).  Plate  85  shows  a  schematic  representa¬ 
tion  of  the  gate  in  horizontal  position.  With  reference  to  the  right-handed  set 
of  axes  XYZ,  the  coordinates  of  the  centroidal  point  of  the  torque  tube  are 
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given  by  CT  (-2.8917,+ 15.6000,+ 1 1.9333),  and  the  coordinates  of  the  cen- 
troidal  point  for  the  gate-torque-tube  system  are  given  by  CG(-1.0390, 

+ 15.200,  + 1 1 .9333).  The  total  mass  of  the  gate-torque-tube  system  is  given 
by  M  =  0.178  Ib-sec^/in.,  whereas  the  tensor  of  inertia  is  given  by  (Boresi 
and  Sidebottom  1985) 


[/] 


76.7 

-2.27 

-2.20 

-2.27 

38.3 

32.2 

-2.20 

32.3 

45.1 

lb-in. -sec^  . 


Applying  twice  the  parallel-axis  theorem  to  transform  (45.1  Ib-in.-sec^) 
into  one  obtains  (Boresi  and  Sidebottom  1985) 


=  4  +  4/(-228.53  in.2) 
=  4.42  lb-in. -sec^ 

=  ^o 


as  the  mass  moment  of  inertia  of  the  gate-torque-tube  system  about  the  torque- 
tube  axis,  Z'. 

If  Young’s  modulus  of  elasticity  for  brass  is  taken  as  £  =  15,500  ksi 
(Flinn  and  Trojan  1990),  and  the  moment  of  inertia  of  both  operating  arms 
about  their  neutral  axis  in  bending  is  calculated  as 


2  • 


12 


recognizing  that  L  =  9.3  in..  Equation  36c  provides  an  estimation  of  the  ftin- 
Hamftntal  frequency  of  vibration  of  the  system  as 

fa)2  =  251,380  sec-2 
.-.  (0  =  501.4  sec'* 
or  /  =  79.8  Hz  , 
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which  compares  well  with  the  finite  element  calculation  for  the  first  mode  of 
vibration  of  the  physical  model  for  the  gate  system  with  a  rigid  torque  tube, 
standing  on  three  twist-hinges  on  a  rigid  foundation  (see  Plate  86.  f  — 

76.3  Hz).  This  structural  parameter  is  expected  to  reflect  a  more  rigid  system 
than  the  actual  physical  model  standing  on  Teflon  bearings  and  load-cell 
foundations  (for  this,  /  »  47  Hz.  See  Plate  83b). 

When  the  system  is  subjected  to  the  steady-state  component  of  the  hydrau¬ 
lic  force  on  the  upstream  side  of  the  gate,  the  general  solution  of  Equa- 
tion  36b  is  the  summation  of  the  complementary  and  the  particular  solutions 
(Naudascher  1991,  Clough  and  Penzien  1975),  i.e.. 


(y4sina)yjr+  ^Bcosw^r)  + 


where 


(0^  =  Wl 


is  the  damped  natural  circular  frequency  of  the  system,  and  A  and  B  are  con¬ 
stants  of  integration  to  be  determined  from  the  initial  conditions. 

Assuming  initial  conditions  for  rest,  the  steady-state  response  of  the  gate 
system  is  given  by  (Clough  and  Penzien  1975) 


(1  -  e^*'“'cos(o^/)  . 

(o2/ 


(37d) 


Evidently,  after  some  time,  the  steady-state  response  tends  asymptotically 


Inserting  the  experimental  values  for  and  /?2  obtained  in  tests  T1087  and 
T1109,  properly  scaled  for  the  physical  model  =  15'^),  and  the  values 
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of  Tj  and  r2  given  by  Equations  26a  and  26b,  properly  scaled  for  the  physical 
model  (X^  =  15'*),  Equation  37d  gives 

(f>  ^  5.5  •  10-^  rad  =  3.5*  ,  (3^0 


which  is  essentially  in  agreement  with  the  experimentally  observed  value,  if 
this  is  taken  as  the  mean  between  the  maximum  and  minimum  readings  during 
the  test. 


Random  Vibration  Model  for  Fluctuations 

The  finite  element  solution  of  the  eigenproblem  for  the  elastic  gate  system 
standing  on  flexible  bearings  on  load-cell  foundations  yields  for  the  first  fre¬ 
quency  of  free  vibration  a  value /j  =  46.91  Hz,  which  compares  extremely 
well  with  the  identified  value  from  the  Twang  test  T0137  on  the  physical 
model  in  dry  air:  /j  =  46.875  Hz  (Plate  83b).  The  value  of /,  =  46.9  Hz  is 
taken  as  the  true  value  of  the  natural  frequency  of  vibration  of  the 
instrumented  gate  system  in  dry  air. 

Letting 

/,=  fir  ,  (38a) 


where 

/p)  =  natural  frequency  of  vibration  of  the  rigid  torque-tube  gate  system 
fij-  =  constant  to  be  determined  experimentally,  i.e.. 


46.9  ^  6S 
79.8  25 


(38b) 


one  has 


(38c) 


Now  when  the  gate  with  water  on  its  upstream  side  is  allowed  to  vibrate 
freely,  an  effective  added  mass  of  water  increases  the  gate  rotational  inertia 
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about  the  torque-tube  axis.  In  other  words,  the  natural  frequency  of  vibration 
for  the  upstream-wet  gate  may  be  expressed  as 


f,  _ 

'  ■  2ir  ^j 


3  El 


where 


K  -  A  ■  4 


(38d) 


(38e) 


and  ^lo  is  a  constant  to  be  determined  experimentally. 

Using  the  experimental  value /,'  =  39.0625  Hz  obtained  in  test  T0139 
(Plate  84b),  the  value  of  is  identified  from  Equation  38c  as 

Denoting  by  />,•  the  pressure  acting  at  the  i***  pressure  cell  location  on  the 
upstream  side  of  Ae  gate  (i  =  1,2,... 9),  assum^  as  uniformly  distributed 
over  its  tributary  area  (/4,),  and  denoting  by  the  distance  from  the  line  of 
action  of  the  force  associated  with  to  the  torque-tube  axis,  the  governing 
equation  of  motion  for  the  gate  subjected  to  the  hydraulic  pressure  field  on  the 
upstream  side  of  the  gate  is  obtained  as  (Clough  and  Penzien  1975) 


K,  +  2^^(2u/;')k, +  (27r/;')2K, 


4' 


9 


E 


(AY/)a(0 


(39a) 


where 


'  <i>  and 


(39b) 


Consider  the  contribution  of  only  to  Equation  39a,  and,  in  particular, 
when  Pi  =  5(t),  i.e.,  a  Dirac  delta  function  applied  at  the  starting  time  (Bendat 
and  Piersol  1993,  Yang  1986): 
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ld^\ 

Pj,.  +  2  ^^(27r/;')  K,,  +  •  (/l,Yy)S  (/)  . 

V"'®  / 


By  Fourier  transforming  Equation  39c,  one  obtains  the  complex  frequency- 
response  function  between  the  pressure  /?,•  and  displacement  V]  of  pressure 
cell  1  in  the  direction  locally  normal  to  the  gate  skin: 


'■(f) 


Therefore,  the  Fourier  transform  of  Vj  is  given  by  (Bendat  and  Piersol 
1993,  Crandall  and  Mark  1963,  Yang  1986) 

V„if)  =  H„if)  •  P,(/)  , 


where  is  the  Fourier  transform  of 

The  Fourier  transform  of  the  total  response  Vj(r)  is  conveniently  obtained 
by  superposition  in  the  frequency  domain  as 


v,{f)  =  Y.  •  W  • 


By  taking  the  inverse  Fourier  transform  of  Equation  42,  the  total  response 
Vjif)  in  the  time  domain  is  obtained  as 


=  E  (AYy)  / 


The  total  elastic  resisting  moment  at  the  root  of  the  operating  arms  is  given 
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or 


Pf  i-i  -00  J 


IF^> 


If  the  pressure  is  assumed  a  zero-mean  stationary  process  (by  extracting 
the  time-average  of  the  record),  then  the  response  Vj,-  and  its  contribution  to 
the  elastic  resisting  moment  will  also  be  zero-mean  stationary  processes, 
and  consequently,  the  following  relations  among  their  corresponding  spectral 
density  functions  hold  (Newland  1984,  Crandall  and  Mark  1963,  Bendat  and 
Piersol  1993,  Hurty  and  Rubinstein  1964); 


By  adding  the  contribution  of  all  p,’s  to  the  total  power,  one  obtains  for  the 
mean  square  of  the  elastic  resisting  moment 


■  / 


(46a) 


in  terms  of  the  half-range  spectral  density  function  (7  (^  (Bendat  and  Piersol 
1993). 

A  program  was  developed  to  run  under  MATLAB  (MATH WORKS  1993) 
for  the  numerical  evaluation  of  Equation  46b.  In  the  particular  case  of  test 
T1080,  the  value  of  the  power  adds  up  to 
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Ml  =  2,204  Gb-ft)- 


which  should  be  compared  with  the  directly  measured  experimental  value, 
i.e., 

Mr  = 

=  90  +  512  =  2,691  Gb-ft)2  , 

where 

=  variance  of  the  random  process  M^j,  given  by  the  area  under  the 
corresponding  spectral  density  function  (Plate  87). 

~  process,  given  by  Equation  37b. 

The  mean  square  of  the  directly  measured  process  differs  from  its  value 
derived  from  the  hydraulic  pressure  field  because  the  directly  measured  pro¬ 
cess  is  really  the  bending  moment  at  the  root  of  the  operating  arms,  which  is 
different  from  the  hydraulic  moment  about  the  torque-tube  axis;  however, 
results  compare  generally  well  (8-percent  error). 


(47) 
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Plate  79.  Lateral  view  of  gate  components,  properties,  and  external  forces 
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Plate  81 . 


Mathematical  model.  Structural  coordinates  and  single  degree  of  freedom  A 


Chapter  6  Engineering  Model 


137 


Moment,  kip-ft  Moment,  kip-ft 
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Plate  85.  Rigid-body  torque-tube  system  in  horizontal  position 


Deformed_plot:  ABAQUS-18-Aug-94  15:0 


Power  Spectral  Density 
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Plate  87.  Spectral  density  of  moment  at  root  of  operating  arms 


7  Conclusions  and 
Recommendations 


A  mathematical  model  to  evaluate  the  total  hydraulic  pressure  field  acting 
on  the  upstream  face  of  a  typical  Montgomery  Point  torque-tube  gate  physical 
model  was  developed  and  experimentally  verified.  This  semiempirical  model 
successfully  predicts  the  upstream  pressure  field  acting  on  most  of  the  gate 
skin,  except  over  deep  locations,  where  the  estimation  is  slightly 
unconservative. 

A  subsequent  experimental  program  on  a  1/15-scale  physical  model  led  to 
a  recommended  design  pressure  line  embracing  both  the  theoretical  prediction 
and  all  test  results  involving  the  regular  operating  conditions  of  the  gate.  This 
design  pressure  line  was  used  as  input  to  conduct  finite  element  analyses  on 
two  separate  structural  idealizations:  one  for  the  physical  model  (made  out  of 
brass  and  standing  on  three  supports)  and  another  for  the  corresponding 
extrapolation  for  the  prototype  (made  out  of  structural  steel  and  standing  on 
five  supports).  The  resulting  reactions  at  the  base  of  the  torque-tube  gate  and 
the  corresponding  force  in  the  actuator  cylinder  were  reported  at  the  time 
(First  Data  Report)  for  the  design  of  support  bearings  and  the  actuator, 
respectively. 

Preserving  the  natural  nappe  on  the  gate  in  service  increases  the  effective 
hydraulic  pressure  on  the  upper  skin  surface  but  probably  not  sufficiently  so 
as  to  justify  the  implementation  of  an  externally  attached  device  to  ventilate 
the  nappe  permanently. 

Two  separate  mathematical  models  were  derived  to  provide  the  designer 
with  information  as  to  the  extreme-value  pressure  field  on  the  gate  that  can  be 
expected  to  occur  during  its  useful  life  with  various  degrees  of  uncertainty, 
based  on  the  information  collected  from  typical  experiments  with  a  fixed  unit 
discharge  of  approximately  50  cfs/ft  (Phase  1  of  the  experimental  project). 
Subsequently,  fiirther  tests  in  Phase  2  of  the  experimental  project  (allowing 
variable  unit  discharges  and  generally  different  pool  conditions  from  those  in 
Phase  1,  and  sometimes  involving  improper  operating  conditions)  provided 
even  more  unfavorable  pressures,  thus  giving  the  controlling  forces  for 
design. 


144 


Chapter  7  Conclusions  and  Recommendations 


The  results  of  the  two  phases  of  the  experimental  project  were  in  agree¬ 
ment  to  confirm  that  the  tests  controlling  the  gate  design  generally  involve  the 
condition  that  all  gates  are  simultaneously  fully  raised  at  some  instant  during 
the  experiment.  In  general,  the  dynamic  fluctuations  of  structural  response 
about  the  steady-state  mean  value  were  only  negligible  to  small  in  magnitude 
during  the  experiments. 

The  gate  eigenproperties  were  identified  by  separate  experimental  tech¬ 
niques,  giving  results  in  close  agreement  with  each  other  and  with  the  finite 
element  eigensolution.  Knowledge  of  the  gate  eigenproperties  is  important  in 
order  to  use  correct  parameters  in  future  modal  analyses  of  the  gate,  to  com¬ 
pare  these  parameters  in  relation  to  the  actual  prototype  gate  as  it  is  con¬ 
structed,  and  to  keep  track  of  structural  deterioration  over  the  life  of  the  gate 
in  service,  due  to  weathering,  corrosion,  local  fatigue  of  component  elements, 
or  simply  due  to  normal  wear  (safety  assessment). 

An  engineering  model  based  on  fundamental  principles  of  structural 
dynamics  was  derived  to  estimate  the  response  of  the  gate  system  subjected  to 
the  upstream  hydraulic  pressure  field.  The  model  involves  two  separate  fac¬ 
ets.  The  first  part  presents  the  deterministic  derivation  of  the  steady-state 
component  of  the  gate  response.  The  second  part  presents  a  random  vibration 
model  for  the  stochastic  description  of  the  fluctuations  of  structural  response 
about  the  steady-state  component.  Both  methods  compare  satisfactorily  with 
the  corresponding  experimental  verification. 

In  an  elemental  spectral  analysis  of  the  test  records,  the  study  evaluated 
the  safety  of  a  typical  torque-tube  gate  in  service  against  fatigue  failure  from 
flow-induced  vibrations  over  the  long  run.  These  initial  studies  indicate  no 
significant  flow-induced  vibration  problems  in  the  model,  since  the  observed 
flow  conditions  produced  no  outstanding  oscillation  energy.  Therefore,  no 
evidence  was  detected  pointing  to  flow-induced  vibrations  of  intensity  such  as 
to  produce  fatigue  in  components  of  the  prototype;  however,  this  does  not 
ensure  categoric  safety  against  fatigue  in  the  prototype. 
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Appendix  A 
Spectral  Analysis 


In  this  appendix,  the  outstanding  issues  resulting  from  the  spectral  analysis 
of  the  pressure  field  data  are  discussed.  Naturally,  only  selected  spectra  for 
the  controlling  experimental  pressure  records  are  actually  exhibited.  Direc¬ 
tions  for  future  studies  to  assess  the  fatigue  performance  of  the  gate  system  in 
the  long  term  are  identified,  but  a  thorough  study  of  this  important  aspect  of 
the  gate  safety  deserves  a  separate  research  project  in  itself 

Tests  T1080  and  T1085,  conducted  at  a  sampling  rate  of  500  sps,  differ 
only  in  the  duration.  Test  T1080  lasts  100  s  whereas  test  T1085  lasts  twice  as 
long.  Plate  A1  shows  the  pressure-cell  locations  for  the  selected  pressure- 
record  spectra  actually  exhibited  in  Plates  A2  through  A7.  These  plates  show 
the  spectral  densities  for  the  upstream,  downstream,  and  net  pressure  records 
in  sequence,  at  each  location  and  for  both  tests.  Each  power  spectral  density 
function  appears  within  the  95-percent  confidence  spectral  band  for  the 
frequency  domain  0  <  f  <  250  Hz.  The  plates  indicate  that  the  test  duration 
has  very  little  effect  on  the  spectral  characteristics  of  the  data. 

Observing  the  upstream  and  the  net  pressure-record  spectra,  one  notices 
that  the  most  important  amounts  of  energy  appear  concentrated  in  the  region 
10  <  f  <  50  Hz  (with  the  main  peak  at  around  40  Hz)  and  in  a  secondary 
energy  pocket  in  a  region  around  60  Hz.  These  frequencies  correspond 
approximately  to  the  first  two  natural  frequencies  of  the  physical  model  with 
upstream  water  only  (see  Plate  83b  for  the  exact  location  of  the  fundamental 
frequency  of  free  vibration).  In  fact,  the  frequency  content  of  the  flow  exhib¬ 
its  peaks  at  20,  30,  40,  and  60  Hz  (Plates  A8  and  A9,  for  test  T1087),  and 
therefore,  clearly  tends  to  excite  the  first  (rigid-plate  rotation  about  the  torque- 
tube  axis)  and  second  (rigid-plate  twist)  modes  of  vibration  (A  similar  phe¬ 
nomenon  has  previously  been  reported  in  the  literature  (March  and  Elder 
1992)).*  The  energy  (area  under  the  spectrum)  associated  with  these  peaks  is 
generally  small-to-moderate  and,  for  the  major  part,  the  flow  spectrum  tends 
to  be  white  noise  in  character,  at  low  energy  levels  (notice  also  that  the  pres¬ 
sure-cell  records  are  influenced  by  the  vibration  of  the  gate  itself,  since  they 
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are  attached  to  the  gate  skin);  nevertheless,  attention  should  be  paid  to  this 
trend  in  a  more  complete  fatigue  analysis  of  the  gate  system. 

Notice  that  there  are  also  energy  pockets  in  the  frequency  regions  around 
0  Hz  and  150  Hz.  The  first  pocket  corresponds  to  the  so-called  DC- 
component  of  the  input  signal  associated  with  the  mean  value  of  the  hydraulic 
pressure  (Bendat  and  Piersol  1993).  The  second  pocket  (~  150  Hz)  tends  to 
excite  higher  modes  of  vibration  and  is  considered  inconsequential,  because 
the  amount  of  energy  associated  with  the  pocket  is  very  small. 

Tests  T1086  and  T1087  were  both  conducted  at  a  sampling  rate  of  100  sps 
for  500  s.  However,  test  T1086  involves  an  integral  nappe,  whereas  in  test 
T1087  the  nappe  is  ventilated.  Plate  A8  shows  the  pressure-cell  locations  for 
the  selected  pressure-record  spectra  actually  exhibited  in  Plates  A9  through 
A14.  The  same  phenomena  discussed  above  appear  in  these  plates.  How¬ 
ever,  the  ventilation  of  the  nappe  in  test  T1087  substantially  dissipates  the  ten¬ 
dency  toward  resonance  in  the  first  two  modes  of  vibration  in  test  T1086, 
where  large  energy  pockets  in  the  region  around  20,  30,  and  40  Hz  are  rein¬ 
forced  by  similar  frequency  contents  in  the  downstream  pressure  field.  These 
energy  concentrations  at  periodic  bands  in  the  frequency  domain  are  believed 
to  be  caused  by  a  cyclic  interaction  among  the  nappe,  the  suctioning  air 
entrapped  by  the  nappe  and  a  turbulent  downstream  water  pool.  These  obser¬ 
vations  are  confirmed  by  tests  T1108  (integral  nappe)  and  Til 09  (ventilat^ 
nappe),  running  both  at  a  sampling  rate  of  500  sps  during  200  s  and  exhibited 
in  Plates  A15  through  A19. 

Plates  A20  through  A23  illustrate  the  coherence  among  selected  pressure¬ 
cell  records  in  test  T1080,  as  an  example.  The  coherence  is  clearly  perfect 
=  1)  in  the  frequency  range  around  the  first  mode  of  vibration  of  the  gate 
(40  Hz)  among  the  pressure  records  for  cells  1,  4,  and  7  and  tends  to  pick  up 
in  the  frequency  bands  around  the  second  and  higher  modes,  indicating  high 
signal-to-noise  ratios  and  that  these  records  are  a  response  to  the  same  excita¬ 
tion  (Bendat  and  Piersol  1993).  This  effect  is  more  clear  among  the  pressure 
records  for  cells  7,  8,  and  9  and  is  consistent  with  the  magnitude  and  phase 
diagrams  for  the  corresponding  cross  spectra. 

This  initial  study  indicates  no  significant  flow-induced  vibration  problems 
in  the  model,  since  the  observed  flow  conditions  produced  no  outstanding 
oscillation  energy.  Therefore,  no  evidence  was  detected  pointing  to  flow- 
induced  vibrations  of  intensity  such  as  to  produce  fatigue  in  components  of  the 
prototype;  however,  this  does  not  ensure  categoric  safety  against  fatigue  in  the 
prototype. 


A2 


Appendix  A  Spectral  Analysis 


Appendix  A  Spectral  Analysis 


A3 


Plate  A1.  Pressure  cells  in  spectral  analysis  of  tests  T1080  and  T1085 


Plate  A2.  Spectral  densities  for  the  records  of  pressure  cell  4  in  tests  T1080  and  T1085 
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Plate  A3.  Spectral  densities  for  the  records  of  pressure  cell  5  in  tests  T1080  and  T1085 
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Plate  A5.  Spectral  densities  for  the  records  of  pressure  cell  7  in  tests  T1080  and  T1085 
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Plate  A6.  Spectral  densities  for  the  records  of  pressure  cell  8  in  tests  T1080  and  T1085 
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Plate  A7.  Spectral  densities  for  the  records  of  pressure  cell  9  in  tests  T1080  and  T1085 
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Plate  A8.  Pressure  cells  in  spectral  analysis  of  tests  T1086  and  T1087 


Plate  A9.  Spectral  densities  for  the  records  of  pressure  cell  4  in  tests  T1086  and  T1087 
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Plate  A10.  Spectral  densities  for  the  records  of  pressure  cell  5  in  tests  T1086  and  T1087 
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Plate  A12.  Spectral  densities  for  the  records  of  pressure  cell  7  in  tests  T1086  and  T1087 
A14  * _ 
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Plate  A13.  Spectral  densities  for  the  records  of  pressure  cell  8  in  tests  T1086  and  T1087 
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Plate  A14.  Spectral  densities  for  the  records  of  pressure  cell  9  in  tests  T1086  and  T1087 
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Plate  A15.  Pressure  cells  in  spectral  analysis  of  tests  T1 108  and  T1 109 


Plate  A17.  Spectral  densities  for  the  records  of  pressure  cell  7  in  tests  T1 108  and  T1 109 
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Plate  A1 8.  Spectral  densities  for  the  records  of  pressure  cell  4  in  tests  T1 108  and  T1 109 
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Plate  A1 9.  Spectral  densities  for  the  records  of  pressure  cell  1  in  tests  T1 108  and  T1 109 
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Plate  A20.  Pressure  cells  in  coherence  analysis  of  test  T1080 


Plate  A22.  Coherence  between  downstream  pressure-cell  records  in  test  T1O8O 
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Plate  A23.  Coherence  between  net  pressure-cell  records  in  test  T1080 
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